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ABSTRACT 

Combining radio observations with optical and infrared color selection - 
demonstrated in our pilot study to be an efficient selection algorithm for finding 
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red quasars - we have obtained optical and infrared spectroscopy for 120 ob- 
jects in a complete sample of 156 candidates from a sky area of 2716 square 
degrees. Consistent with our initial results, we find our selection criteria - 
J — K > 1.7, R — K > 4.0 - yield a ~ 50% success rate for discovering quasars 
substantially redder than those found in optical surveys. Comparison with UVX- 
and optical color-selected samples shows that > 10% of the quasars are missed 
in a magnitude-limited survey. Simultaneous two-frequency radio observations 
for part of the sample indicate that a synchrotron continuum component is ruled 
out as a significant contributor to reddening the quasars' spectra. We go on 
to estimate extinctions for our objects assuming their red colors are caused by 
dust. Continuum fits and Balmer decrements suggest E(B — V) values ranging 
from near zero to 2.5 magnitudes. Correcting the i^-band magnitudes for these 
extinctions, we find that for K < 14.0, red quasars make up between 25% and 
60% of the underlying quasar population; owing to the incompleteness of the 
2MASS survey at fainter i^-band magnitudes, we can only set a lower limit to 
the radio-detected red quasar population of > 20 — 30%. 

Subject headings: dust, extinction; quasars: general; surveys 

1. Introduction 

The definition and physical understanding of quasars has evolved and expanded since 
the first quasars were found over 40 years ago. Discovered in the radio, quasars were thought 
to be radio-loud objects until Sandage (1965) uncovered a population of objects with optical 
properties identical to quasars, but without radio emission. It is now known that there are 
roughly ten times more radio-quiet than radio-loud quasars. 

The spectroscopic properties of the earliest known luminous quasars showed blue optical- 
to-ultraviolet colors and broad emission lines, as well as narrow forbidden lines of highly 
ionized species. This, coupled with their compact morphologies and extreme luminosities, led 
to the commonly accepted model for the quasar as powered by accretion onto a supermassive 
black hole. This picture was expanded by invoking orientation-based arguments to include 
narrow emission-line spectra of AGN (e.g., Seyfert2 and narrow-line radio galaxies) as well 
as highly- variable featureless spectra (e.g., BL Lacs). 

Such orientation models depict an axisymmetric geometry with an accretion disk at the 
center surrounded by broad emission-line clouds. Dusty, molecular clouds, often modeled 
as a torus coplanar with the accretion disk, obscure the broad emission lines along certain 
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lines of sight. Beyond the obscuring clouds lies the narrow emission-line region (Urry & 
Padovani 1995). This model is successful at explaining the different varieties of AGN spectra: 
whether broad lines are seen depends on our viewing angle toward the AGN. It has also been 
successful at predicting AGN number counts in the infrared as well as explaining the shape 
of the extragalactic x-ray background (e.g., Madau et al. 1994; Comastri et al. 1995; Gilli 
et al. 2001; Ueda et al. 2003; Treister et al. 2004; Treister & Urry 2005; Treister et al. 2006). 

Although the unification scheme is well-established for the low-luminosity AGN, it has 
been problematic for quasars. Several individual cases of Type 2 quasars at high redshift 
(z > 3) were found in the X-ray. These sources are found to have high column densities 
(Nh > 10 23 cm -2 ) and X-ray luminosities (Lx > 10 44 erg s" 1 ) as well as narrow permitted 
emission-lines in their optical spectra and extremely weak optical continua (Norman et al. 
2002; Stern et al. 2002). Zakamska et al. (2003) have presented low-redshift Type 2 quasar 
population by taking advantage of the large spectral database in the SDSS. Selecting for 
objects with 0.3 < z < 0.83, Zakamska et al. (2003) estimate that 50% of their sample (~ 150 
out of 291 Type 2 AGN) have quasar luminosities (Mb < —23). Their multiwavelength and 
polarization properties also support their Type 2 classification and the unification scheme 
for luminous AGN out to the quasar regime (Zakamska et al. 2004, 2005). 

Several observational findings linking supermassive black holes and quasar activity to 
their host galaxies, however, are not accounted for by these orientation models and may 
provide clues about joint formation and evolution of quasars and their hosts. Most galaxies 
harbor a quiescent supermassive black hole, suggesting that a quasar may be a short-lived 
phase in the life of a galaxy. Moreover, the black hole mass appears to correlate with the host 
galaxy's bulge mass and stellar velocity dispersion (Magorrian et al. 1998; Ferrarese & Merritt 
2000). This, combined with the remarkably similar evolution of quasar activity and star 
formation over the history of the universe, suggests that the quasar is an essential ingredient 
in galaxy formation. Initial results from semianalytical models successfully reproduce some 
of these observational trends (Kauffmann & Haehnelt 2000) and more recent numerical 
simulations of joint black hole/galaxy formation have been able to account for many of these 
observations (Di Matteo et al. 2005; Hopkins et al. 2006). In addition, these models predict 
an obscured phase for the nascent quasar as a galaxy merger triggers star formation and 
funnels gas into the nucleus. Eventually, feedback from the quasar blows the dust and gas 
away, quenching the accretion and ending the black hole's active phase (see also Sanders 
et al. 1988; Yun et al. 2004). 

Although reddened quasars are predicted in both the static, orientation-based model 
and the evolutionary scenarios, this population has only begun to be uncovered. The major- 
ity of quasars described in the literature have been found from optically selected samples and 
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have remarkably homogeneous spectra with blue optical-to-ultraviolet spectral energy dis- 
tributions. The recent discovery of a population of red quasars has challenged the view that 
all quasars are blue, UV-luminous sources and suggests that a large population of quasars 
with red colors exists (Webster et al. 1995; Cutri et al. 2001; Gregg et al. 2002; Richards 
et al. 2003; White et al. 2003b; Glikman et al. 2004). The advent of large near-infrared 
detectors that perform at the level of optical CCDs has enabled surveys for quasars that 
are largely insensitive to dust extinction. Since most near-infrared sources in the sky are 
Galactic stars, radio surveys offer an effective complement to near-infrared quasar searches 
by selecting predominantly extragalactic sources, as well as being completely insensitive to 
dust extinction. 

In Glikman et al. (2004, hereafter, Paper I), we used this principle to search for red 
quasars by studying radio sources that are bright at 2 /iin (K < 16) but undetectable at 
R > 20. Spectroscopic followup of these sources revealed that objects with J — K > 1.7 
and R — K > 4 had a 50% chance of being highly reddened quasars. Seventeen red quasars 
were detected in that survey, and we concluded that red quasars make up > 20% of the total 
quasar population. 

In this paper, we expand our pilot survey to produce a complete catalog of objects 
obeying our optical-to-near-infrared color criteria over the same 2716 deg 2 overlap region 
between the VLA FIRST and 2MASS surveys 1 . We review the input catalogs and the 
candidate selection process in Section 2. We describe the spectroscopic followup of our 
candidate list as well as the classification of the spectra and the resultant catalog statistics 
in Section 3. In Section 4, we present the observed surface density of red quasars on the sky 
and compare it to that for optically selected quasars. In Section 5, we derive and compare 
extinctions from continuum fits to a reddened quasar template (§5.1) and extinctions derived 
from Balmer decrements (§5.2). In Section 6 we examine the radio spectral indices of the 
red quasar population. We estimate the percentage of red quasars in the overall quasar 
population in Section 7, while Section 8 discusses variable blue quasars found in our survey. 
Section 9 places our sample in context with other radio-selected reddened quasars in the 
literature. Section 10 summarizes our findings. 

2. The Selection Process 

Our red quasar survey relies upon three wavelength regions to construct our candidate 
list. We require that our quasars be radio sources to avoid confusion with Galactic stars. 



In this work we use the 2000 July FIRST catalog to maintain continuity with Paper I 
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Additionally, since we are looking for red objects, we require that our quasars be bright at 
2 /im. Finally, we use optical data to select objects with red optical-to-near-infrared colors. 

The FIRST Survey (Becker et al. 1995) has mapped ~ 9000 deg 2 of the sky at 20 cm 
using the Very Large Array (VLA) in the B-configuration. Coverage includes ~ 8400 deg 2 
in the north Galactic cap and ~ 600 deg 2 in the south Galactic cap. With the array in this 
configuration, the survey is roughly analogous to all-sky optical surveys, with 5" resolution 
and 0Y5 positional accuracy. The 3-minute snapshot integration time yields a typical rms of 
0.15 mJy The July 2000 catalog generated from this survey has ~ 722, 354 sources brighter 
than the survey's detection threshold of 1 mJy over 7988 deg 2 . 

The 2MASS survey (Skrutskie et al. 2006) used two telescopes, one at Mt. Hopkins, 
Arizona and the other at Cerro Tololo, Chile, to image the sky in three near-infrared bands, 
J (1.24 /xm), H (1.66 /im), and K s (2.16 /mi), between 1997 and 2001. This four-year effort 
has produced a Point Source Catalog (PSC) with ~ 4.70 x 10 8 sources. The survey is 99% 
complete to a magnitude limit at 10a of J = 15.8, H — 15.1 K — 14.3. Fainter sources are 
included in the catalog (down to K < 16) but with lower signal-to-noise and completeness. 
The survey has a resolution of 4" and an astrometric accuracy better than / .'5. We utilize 
the 2nd Incremental Release of 2MASS, which includes observations through 1999, covers 
19,681 deg 2 , and contains > 1.6 x 10 8 point sources 2 . 

The Guide Star Catalog II (GSC-II) is an all-sky optical catalog produced by scanning 
the second-generation Palomar Observatory (POSS-II) and UK Schmidt Sky Surveys at 1" 
resolution and generating a catalog with positions, magnitudes and morphological classifi- 
cations. We employed the GSC 2.2.1 which is a pre-release version not available outside the 
STScI that reaches the POSS-II plate limits of F < 20.80 and J < 22.50. It is similar to the 
publicly released catalog, version 2.3.2. The astrometric accuracy of GSC-II is better than 
1". 



2.1. The FIRST-2MASS Sample 

We began our search for a population of optically obscured, red quasars by looking for 
objects detected in the FIRST and 2MASS catalogs but which lacked an optical counterpart 
in the Cambridge Automated Plate Measuring Machine (APM) scans of the first-generation 



2 Although the 2MASS All-Sky Catalog is described as being of superior quality by the 2MASS Explana- 
tory Supplement (http://www.ipac.caltech.edu/2mass/releases/allsky/doc/explsup.html, §I.6.v) , we discuss 
the subsample selected from the 2nd Incremental Release Point Source Catalog in this work for continuity 
with Paper I and because of our high level of spectroscopic completeness over this well-defined area. 
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Palomar Observatory Sky Survey (POSS-I) plates. Spectroscopic followup of 50 of the 
resulting 69 candidates allowed us to define a region in color-color space which maximized 
our efficiency for finding red quasars. We found that ~ 50% of the spectroscopically identified 
objects in our initial sample defined by R — K 3 > 4 and J — K > 1.7 were heavily obscured 
quasars. The other ~ 50% were galaxies exhibiting some form of activity such as a starburst 
or narrow-line AGN activity. Objects with bluer J — K colors tended to be late-type stars 
and objects with R — K < 4 tended to be elliptical galaxies; these results were reported in 
Paper I. 

Based on these findings we chose to impose these color cuts on FIRST-2MASS sources, 
employing optical data from GSC-II to generate a sample of candidate red quasars. As was 
shown in Paper I, a 2" search radius between FIRST and 2MASS yields a 94% probability 
that the two sources are physically associated. A match between FIRST and 2MASS using 
this search radius yielded 12,435 matches, keeping only the closest match (a catalog we refer 
to henceforth as F2M). We matched these sources, using their FIRST positions, to the GSC- 
II catalog with a search radius of 2". There were 10,903 matches within 2" of the FIRST 
positions in F2M. 

In order to apply the optical color cut to the optical matches we need to take into 
account sources that have been detected in one band but not in another. We replaced 
the magnitudes for undetected objects with their respective plate limits, B = 22.50 and 
R = 20.80. We then applied our J — K > 1.7 and R — K > 4 color cuts, and reduced the size 
of our candidate list to 146 sources. In cases where bright objects were detected on plates 
observed through the Photographic J bandpass but were undetected in Photographic F (e.g., 
saturated stars), artificially red colors may have been assigned. We examined finding charts 
of all 35 sources with R > 20.80 and found that 5 were indeed too faint to be discerned 
on the Digitized Sky Survey images; we kept them in the candidate list. We obtained 
magnitudes for the remaining 30 sources using auxiliary catalogs. Twenty-six objects had 
detections in the SDSS-DR5 catalog which contributed r magnitudes - all but three were 
too blue to remain in our candidate list. We obtained APM E magnitudes for the remaining 
4 sources and removed all but one object, F2M J171058. 4+310958, which has E = 19.34 and 
E-K = 4.38. The final list of FIRST-2MASS sources with GSC II matches within 2" thus 
contains 120 objects with R — K > 4 and J - K > 1.7. 

In addition, fourteen objects too faint to be detected in the GSC-II catalog were previ- 
ously reported in Paper I as having J — K > 1.7. We include these sources in our candidate 



3 B and R are used to refer to "blue" and "red" filters in cither the POSS-I E and O or the POSS-II F 
and J bands as well as the SDSS r band. 
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list, adopting the plate limits as lower limits to their B and R magnitudes. 

The 2MASS survey quotes upper limits for objects with a photometric signal-to-noise 
ratio smaller than ~ 3. Objects that are below the 2MASS J-band limit can have J — K 
values that are too blue to be selected into our sample, though they might have J — K > 1.7. 
There are 718 F2M sources with a J-band signal-to-noise ratio less than 3 and J — K < 1.7. 
618 of these have a GSC-II match within 2" and 90 have no GSC-II match . In addition, we 
checked the 2MASS All-Sky catalog whose photometry is based on a different thresholding 
algorithm and found four J-band detections. One of these had J-K=1.52, and we remove it 
from the candidate list. We examined the POSS-II images of the optically unmatched objects 
and removed 72 large galaxies and saturated stars which made it into this list as a result of 
bad photometry in at least one survey. We applied the R — K color cut to the objects with 
GSC-II matches and were left with 20 sources. We obtained SDSS r magnitudes for 15 of the 
18 optically unmatched sources and APM photometry for the other 3. Two sources survive 
the R — K color cut. Therefore, in order for our survey to be a i^-band limited survey, 
rather than a J-band limited survey (which would be more susceptible to dust extinction) 
we append to our candidate list 22 objects with a J-band signal-to-noise ratio smaller than 
3 that have nominal J — K values less than 1.7, but have limits consistent with J — K > 1.7. 
This may contaminate our candidate sample slightly, but assures completeness. 

The final candidate list thus contains 156 sources, 140 of which have detections in the 
GSC-II catalog, and 16 of which are fainter than the sensitivity limit of GSC-II. This list 
includes the 31 objects from Paper I which had no detection in the APM catalog and meet 
our color criteria. We present these candidates in Tables 1 and 2 (which we describe in §3.2). 
Table 1 lists the sources that have J-band detections in 2MASS and Table 2 lists the 2MASS 
objects whose J-band detection had a signal-to-noise ratio less than 3 and J — K < 1.7. 

The 2nd Incremental Release of 2MASS is made up of a noncontiguous set of tiles over 
the celestial sphere, while the FIRST covers a solid "footprint" over the North and South 
Galactic caps. We computed the area of the F2M survey in Paper I by counting the number 
of FIRST sources that fell within the four corners of a 2MASS tile. Since FIRST sources 
are smoothly distributed over the sky, the fraction of FIRST sources inside the 2MASS area 
determined the fraction of the FIRST area covered by the two surveys. This amounted to 
34%, or 2716 deg 2 . 

Figure 1 shows a schematic diagram of our selection process. The boxes with dotted-line 
borders represent the input catalogs: FIRST, 2MASS and GSC II. The matching radius is 
given in a circle at the node of each match. The final catalog is indicated by four bold- 
lined boxes, one representing the 120 sources with FIRST, 2MASS and GSC II matches and 
the colors R - K > 4 and J - K > 1.7, a second for the 14 FIRST-2MASS sources with 
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no GSC-II detection, and two final boxes indicating the 22 additional sources with J-band 
upper limits. 

3. Spectroscopic Followup 
3.1. Observations 

Spectroscopic followup of our candidates was carried out over a five-year period at 
four observatories at both optical and near- infrared wavelengths. A majority of the optical 
spectra were obtained at the 10 m W. M. Keck Observatory in Hawaii using the ESI and, to 
a lesser extent, the LRIS spectrographs. A handful of spectra were obtained with the Kast 
spectrograph on the 3 m Shane telescope at Lick observatory. We also obtained spectra from 
publicly available databases such as the FBQS II and III (White et al. 2000; Becker et al. 
2001) as well as the Sloan Digital Sky Survey Spectral database. The LRIS spectra were 
reduced using standard IRAF procedures, while the ESI spectra were reduced using both 
IRAF tasks and software designed by the authors (see White et al. 2003a). 

The near-infrared observations, spanning 0.8 — 2.5 /im, were conducted at the NASA 
Infrared Telescope Facility (IRTF) with the SpeX spectrograph (Rayner et al. 2003). We 
reduced these spectra using the Spextool software which is specially designed for SpeX data 
(Cushing et al. 2004). We corrected for telluric absorption using the spectrum of a nearby 
A0V star obtained immediately before or after each object spectrum using the technique 
outlined in Vacca et al. (2003). 

Five near-infrared spectra, spanning 0.95 — 2.5 /im, were obtained at the MDM Ob- 
servatory on Kitt Peak with the TIFKAM infrared camera in spectroscopic mode. These 
spectra were reduced using standard IRAF procedures. We corrected these spectra for tel- 
luric absorption using the spectrum of a nearby AO V star as well as a G2 V star following 
the technique outlined in Hanson et al. (1996). 

3.2. Results 

We have obtained spectra for 120 of our 156 candidates from optical and/or near-infrared 
observations or from data found in the literature. Our complete candidate lists are presented 
in Tables 1 (for objects detected in the 2MASS J-band) and 2 (for objects with J — K < 1.7 
but with a J-band signal-to-noise ratio less than 3 which, as a result, may result in a true 
vale of J - K > 1.7). The FIRST catalog J2000.0 Right Ascension and Declination is listed 
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in columns (1) and (2) followed by the FIRST peak and integrated 20cm flux densities in 
columns (3) and (4). We list the GSC-II B and R A magnitudes in columns (5) and (6) 
and the 2MASS J, H, and K s magnitudes in columns (7), (8), and (9), respectively. We 
list the J — K and R — K colors used in the selection of this sample (columns 10 and 
11). When available, we present the object's redshift and classification (described below) in 
columns (12) and (13). In column (14) we list a reference for the optical and/or near-infrared 
spectrum used to classify our source. Column (15) provides any references to the object in 
the literature. 



3.2.1. Object Classification and Narrow-Line AGN 

We classified our spectra using continuum and line analysis. We cross-correlated our 
spectra with a set of galaxy templates from Kinney et al. (1996). We confirmed the classifi- 
cation of spectra with no obvious emission lines that were well-correlated with an elliptical 
or SO template by measuring the 4000A break at the same redshift. In spectra with emission 
lines, we measured the redshifts and line profiles, classifying objects with line widths > 1000 
km s" 1 as quasars. We use the redshifts determined from the line-measurements to confirm 
the redshifts determined from the cross correlations. 

Twenty-five spectra had narrow, permitted emission lines. We classified twenty of them 
by fitting Gaussian profiles to the lines and performing line diagnostics following Kewley 
et al. (2006). Figure 2 and Table 3 list the results of these measurements. Seven objects 
were classified as emission-line galaxies and ten objects were classified as AGN in all three 
diagrams. We assigned the classification that agreed in two diagrams for the remaining 
sources. The redshifts of two more objects were too high to measure Ha. Their [O IIl]/H/3 
ratios, however, are ~ 1 making them very likely AGN. The three remaining narrow-line 
spectra have only a near-infrared spectrum revealing narrow Pact and Pa/3. We label these 
as narrow-line emitters ("NLemit") in Table 1. 

The classification of the 120 spectra breaks down as follows: Thirty-seven galaxies 
(ten starburst galaxies and twenty-seven red elliptical galaxies); fourteen narrow-line AGN 
(including one Seyfert 1.5 from Smith et al. 2002) and three narrow-line emitting spectra 
that could not be further classified; one star; five BL Lacs and 57 quasars (line widths > 
1000 km s _1 ) of which 52 are red quasars and 5 are blue quasars (discussed in more detail 
in §8). We also obtained two IRTF spectra and two optical spectra which we were unable 
to classify due to poor signal-to-noise and/or a lack of identifiable features. We assume that 



4 Photographic F and J magnitudes. 
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these objects are not quasars in our subsequent analyses (although some of them could be). 
Our efficiency for finding broad-line quasars holds at ~ 50% (57/120) as we predicted in 
Paper I. 

In this paper, we include only the (type I) quasars in our analysis. Rest-frame visual 
extinction in narrow-line (type II) AGN can reach hundreds of magnitudes and the de- 
reddening techniques that we use (in §5) cannot be applied. Without a reddening estimate, 
we cannot compute the intrinsic luminosity of these sources for comparison with their unred- 
dened counterparts. Furthermore, the nuclear light in type II AGN is entirely obscured even 
in the near-infrared. The near-infrared colors of type II AGN are consistent with coming 
from the host galaxy (Zakamska et al. 2004). Since our color-selection technique is designed 
to find lightly reddened quasars, and intentionally avoids galaxies (see Figure 4 of Paper I), 
our sampling of type II AGN is likely very incomplete. Similar distinctions have been made 
by Lacy et al. (2007) and Brown et al. (2006) who estimate the fraction of reddened type I 
quasars separately from type II sources. 

We note that as in Paper I several of our red quasars had all narrow-lines in their optical 
spectra. It was the near-infrared spectrum that revealed broad Pa/3 or Ha. Without the near- 
infrared spectrum, these sources would have been classified as narrow-line AGN. Fourteen 
of the sixteen objects classified as "NLAGN" or "NLemit" have only optical spectra. It is 
possible that near-infrared spectroscopy would reveal broad emission lines for some of these 
objects. Therefore, for this and other reasons discussed below, our estimate for the fraction 
of missed red quasars is only a lower limit. 

We lack a spectrum for one of our quasars, which we identified from the literature. 
This is a ~ 400 mJy source and overlaps with the Webster et al. (1995) red quasar sample 
(Drinkwater et al. 1997). The remaining 56 quasar spectra are presented in Figure 3. We 
have both an optical and near-infrared spectrum for 33 of our quasars. Since the overlapping 
regions between the optical and near-infrared spectra (8000 — llOOOA) were often the noisiest 
portions of each spectrum, the size of the overlap region used in a given pair of spectra covered 
a large range (100A to over 2500A). In addition, the optical and near-infrared observations 
were not contemporaneous and therefore the spliced spectra are susceptible to variability 
effects. We chose to scale the spectra individually using our own interactive software. We 
overplotted the optical and near-infrared spectra and selected a region for each pair in the 
least noisy regions of overlap. We then scaled the spectrum with the lower median flux in 
that region to that with the higher median. Since our data were taken on different nights, 
seeing and transparency variations as well as slit losses resulting from different aperture sizes 
required scaling factors up to 12; most of the spectra (23/33) required scaling by less than 
a factor of 3. 
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We note that we are 100% spectroscopically complete to K < 14.5; our completeness 
drops somewhat to 77% for the entire sample, as a result of our ~ 70% completeness for 
K > 14.5. Figure 4 shows a histogram of the K magnitudes binned by 0.25 magnitudes. 
The unshaded histogram contains all 156 candidates, while the shaded histogram shows all 
identified objects (including the four unclassified spectra, which we include here as non- 
quasars). We overplot the 57 quasars in black. 

4. The Observed Surface Density of Quasars on the Sky 

In order to conduct a meaningful comparison between our red quasar sample and the 
blue quasars found in optical surveys, we draw upon the FIRST Bright Quasar Survey 
(FBQS; Gregg et al. 1996). The FBQS looked for quasars in matches between FIRST radio 
sources and sources in the APM scans of the first generation Palomar sky survey (POSS-I) 
plates (separations < l'/2) which were classified as point sources on either the photographic E 
(red) or O (blue) plates. The initial catalog contains 636 quasars over 2682 deg 2 and imposed 
an E = 17.8 magnitude limit on the APM sources as well as an O — E < 2 color criterion 
(FBQS II; White et al. 2000). A deeper segment of the survey covered 589 deg 2 and found 
321 quasars with the same color cut, but with a deeper magnitude limit of E = 18.9 (FBQS 
III; Becker et al. 2001). Since we are interested in objects detected in both FIRST and 
2MASS, we matched the FBQS II and III catalogs to the 2MASS All-Sky Data Release with 
a search radius of 2" and found 503 and 134 matches to FBQS II and III, respectively. These 
quasars form an optically selected radio sample in the same .fT-band flux-limited survey and 
can be compared to our sample. 

In Figure 5 we plot the B—K color distribution of FBQS II and III quasars along with the 
F2M red quasars as was done in Paper I. This sample now fills in the 5 < B — K < 6.5 region 
of the plot which was previously unpopulated since our earlier optical selection criterion 
required objects to have B > 21.5. 

The redshift distribution is shown in Figure 6. Our survey is able to reach a redshift 
range similar to the FBQS II and III samples, although the density of our sources decreases 
with redshift. This is a consequence of the increased effect of dust on observed if-band light 
as a function of redshift. At higher redshifts, the rest-frame optical light is shifted into the 
.fT-band. This light is more susceptible to extinction resulting in an observed magnitude 
that falls below the 2MASS detection threshold for all but the most luminous objects. 

We present the spatial density of red quasars compared with FBQS II, III and a UVX- 
selected sample in Table 5. The UVX-selected sample is described in Paper I. Briefly, it 
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was constructed by combining the Large Bright Quasar Survey (LBQS) described in Hewett 
et al. (1995) and the Two-Degree Field QSO Redshift Survey (2QZ) described in Croom et al. 
(2001). We matched these sources to FIRST and 2MASS and created a UVX-selected sample 
that is comparable to the F2M red quasars. This sample goes deeper in the optical than the 
FBQS samples. The magnitude range for the LBQS and 2QZ quasars is 16.0 < Bj < 20.85. 

We plot the densities as a function of K magnitude in Figure 7. Assuming our selection 
method finds all obscured quasars bright enough to be detected in the 2MASS and FIRST 
surveys (which is unlikely - see discussion below), the total surface density of quasars ought 
simply to be the sum of the blue quasars plus the red quasars, normalized by their respective 
coverage areas. Depending on the optical survey used for comparison, the last line of Table 5 
reveals that FBQS II misses 12.0±2.1%, FBQS III misses 12.4±2.1%, and the UVX selection 
method misses 8.2±1.4% of all FIRST-detected quasars visible in the 2MASS survey. These, 
of course, are lower limits to the fraction of red quasars at a given magnitude. In addition to 
our spectroscopic incompleteness, the observed if-band magnitudes of our quasars are still 
affected by dust and we have compared counts as a function of unreddened magnitudes with 
magnitudes that need to be corrected for absorption. In the following section we determine 
the reddening in each quasar using its spectrum. We can then correct the observed .fT-band 
magnitudes for absorption and compare their counts to the unobscured samples. 

We note that although we have made an effort to compare samples with as similar 
properties as possible by comparing the F2M red quasars to objects detected in both FIRST 
and 2MASS, the fraction of missed quasars that we compute is a lower limit to the real 
fraction of missed quasars. This is because there are quasars with FIRST and 2MASS 
detections that may not have been selected by any of these surveys. The FBQS II survey 
has an APM E = 17.8 magnitude limit. Any object that is fainter than K = 13.8 (which 
is true of ~ 80% of the F2M quasars) and that is a slightly bluer than R — K = 4 will not 
be in either the FBQS or the F2M red quasar sample. In the FBQS III case (which has an 
E = 18.9 magnitude limit) this is a smaller effect - only quasars with K > 14.9 (amounting 
to about 25% of the F2M quasars) can be missed. Overall, however, this is likely to be a 
small effect, since the space densities of FBQS II and III quasars with 2MASS matches are 
nearly identical (see Figure 7), meaning that the fainter optical magnitude limit of FBQS 
III does not recover significant numbers of slightly reddened quasars missed by FBQS II. 

Taking into account the possible existence of additional missed quasars, the actual 
observed fraction of red quasars is no longer simply the number F2M red quasars divided by 
sum of the blue quasars plus the red quasars, normalized by their respective coverage areas. 
The quasars missed by all these surveys enter in both the numerator and denominator of 
this equation and as a result our computed missed fraction above must be a lower limit to 
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the true missed fraction of quasars. 

Higher fractions of missed reddened type I quasars have been found by Brown et al. 
(2006) and Lacy et al. (2007) by selecting quasar candidates in the dust-insensitive mid- 
infrared. Brown et al. (2006) find that reddened type I quasars make up ~ 20% of the type 
I quasar population, while Lacy et al. (2007) find an even higher fraction of 30%, based on 
a sample of 11 reddened and 25 unreddened type I quasars. When we take reddening into 
account for our sample, we also find a higher missed fraction (see §6). 



5. Reddening Investigation 

5.1. Continuum Fits to a Reddened Quasar Template 

In Paper I we assumed that the red colors of our quasars were due to extinction by dust 
and derived the color excess, E(B — V), in the rest frame for each source. Here, we fit to 
each object a reddened quasar composite from the FIRST Bright Quasar Survey (Brotherton 
et al. 2001) combined with a near-infrared quasar composite (Glikman et al. 2006), resulting 
in wavelength coverage from 800 to 35,000 A. We linearly fit the reddening law parameterized 
by Fitzpatrick (1999), k(\), to the reddening of each individual quasar, /(A), with respect 
to the unreddened composite, /o(A). Assuming standard extinction, 

/(A) = / (A)e-*\ (1) 

our linear fit has the form: 

\nJ f ^\ r k(X)E(B - V) 

log Km) = = — ^08^- (2) 

We exclude strong emission lines and use a robust linear fit which ignores extreme outliers 
in very noisy spectra and returns E(B — V) as the slope. 

Twenty-two of our 55 quasar spectra have only optical or near-infrared spectra and we 
adopt the amount of reddening that gave the best fit for each of these spectra. Figure 8 
shows examples of these fits for two objects, one with an optical spectrum and one with a 
near-infrared spectrum. 

This method of determining E(B — V) does have several drawbacks for determining the 
reddening in objects with combined optical and near-infrared spectral coverage. Our spectra 
were not obtained simultaneously; the interval between observations ranges from several 
weeks to a few years. Aside from variability in brightness, which involves a simple scaling 
of one spectrum to the other, variability in the spectral shape is also likely and cannot be 
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accounted for. This is an especially important caveat in the case of dust-obscured quasars, as 
the obscuring material close to the nucleus may cross the line-of-sight over these timescales 
leading to different measurements of the source's extinction (cf. Risaliti et al. 2005). 

In addition, our choice of scaling described in §3.2 may have a strong effect on the fit. 
Since rest frame UV-optical light is more sensitive to extinction than the near-infrared, we 
expect that the optical spectrum, when available, will provide a more reliable measure of 
the extinction. To determine the effect on our reddening measurements, we fit the reddening 
both to the combined, as well as to the individual optical spectra alone. In most cases, 
fitting to the combined spectra gave satisfactory fits. Figure 9 shows an example of two 
optical-to-near-infrared spectra, and the fitted composite. However, in about a third of the 
cases, the combined fit was extremely poor, while fitting the optical spectrum alone yielded 
an excellent fit. Figure 10 shows an example of such a case. 

Another caveat to fitting a reddened quasar template to our spectra is the presence of 
galactic starlight. Although this is not an issue for luminous, unreddened AGN, where the 
nuclear light outshines the galaxy, in our highly reddened objects the presence of a galactic 
component can have an effect on our fits. This is especially a concern at low redshifts where 
we are sensitive to both intrinsically lower luminosity quasars, and more heavily obscured 
objects. Spectropolarimetry of near-infrared selected quasars from 2MASS confirmed that 
host galaxy light can contribute significantly to the spectra of type II AGN, diluting the 
polarized flux from the AGN (Smith et al. 2003). 

We conducted a two-component fit to our spectra, using the same quasar template plus 
an elliptical galaxy template from Mannucci et al. (2001) allowing the reddening and fraction 
of quasar and galaxy light to be fit simultaneously. The fitting function has the form: 

/(A) = A/ B a(A) + S/ Q so(A) e -^. (3) 

Since the strongest feature in an elliptical galaxy template is the 4000A break, which 
moves out of the optical at z ~ 1.5, and the scaling between the optical and near-infrared 
spectra could introduce additional uncertainty, we used only the optical spectrum in these 
fits. The fits converged and significant galactic components were found only for eight objects 
with z < 0.5. We subtracted the galactic component from the optical spectrum, and re-scaled 
it to the near-infrared spectrum. We then ran the linear fitting procedure (equation [2]) using 
only the quasar template on these objects. We compared E(B — V) derived from the three 
techniques: (1) the original linear fits to the combined spectrum, (2) the two-component 
fits to the optical spectrum and (3) the galaxy-subtracted, linear fits to the recombined 
spectrum. The results generally found the highest derived extinctions using method (2) 
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followed by method (3). 

The following examples illustrate these findings. For F2M170802.4+222725, at z = 
0.377, we derive E(B — V) — 1.05 using technique (1), E(B — V) = 1.57 using technique 
(2), and E(B - V) = 1.10 using technique (3). F2M130700.6+233805 has only an optical 
spectrum, so scaling is not a concern. We derive E(B — V) = 0.58 using technique (1) for this 
object. Technique (2) measures E(B — V) = 0.55 and technique (3) yields E(B — V)= 0.56. 

These results underscore the difficulty of estimating reddening in quasars spanning such 
a large redshift range, whose data are highly heterogeneous. The nature of our sample and of 
the data introduces unconstrained degrees of freedom into our fitting techniques, including 
the optical to near-infrared spectrum scaling, possibly variable reddening as well as other 
variability, and an unknown fraction of galactic starlight contamination. 

Taking all of this into account, while striving for consistency, we disregard the galactic 
contribution and adopted the amount of reddening from either the combined optical and near- 
infrared spectra or just the optical spectrum. Whichever gave the best fit upon examination. 
The measured extinctions range from E(B — V) — —0.05 to 1.5. We deredden the spectra 
to compute their intrinsic if-band magnitudes. 

Table 4 presents the reddening parameters derived for our spectra. We list each object's 
J2000 coordinates in columns (1) and (2) followed by the redshift in column (3) and the 
.fT-band magnitude in column (4). Column (5) lists the i^-band extinction based on the 
E(B — V) value reported in column (9). Column (6) lists the object's R — K color. We 
compare reddenings derived from two different methods in the next three columns. Columns 
(7) and (8) list the E(B — V) values derived using Balmer decrements in the total hydrogen 
line profile and only the broad component of the line profile, respectively, for eligible spectra 
(see §5.2), while column (9) lists the reddening from our composite fitting technique described 
in this section. We list in column (10) the spectrum used to measure the extinction. 

5.2. Reddening Determined from Balmer Decrements 

The line ratios seen in our spectra can provide useful information about the reddening 
processes along the line of sight. Many of our objects have been sufficiently extinguished such 
that the broad emission lines are no longer present in the rest frame optical; it is only from 
broad Paschen lines that we are able to classify them as quasars. In some cases the Balmer 
line profiles show a weakened broad component with sharp narrow lines superimposed. In 
these cases we can measure the extinction in the broad and narrow components of the Balmer 
lines separately to determine if the dust is located behind or in front of the narrow line region. 
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We can also correlate the extinctions measured directly from Balmer decrements with the 
extinctions derived from our continuum fits. 

Forty-five quasars in our sample have spectra with both Ha and H/3 lines included in our 
wavelength coverage. Of these, nineteen have both lines in a single optical or near-infrared 
spectrum. To avoid uncertainty in the Balmer decrements from the scaling of the infrared 
spectrum to the optical, we measure Balmer decrements for only these nineteen quasars. We 
follow the line-fitting technique outlined in Greene & Ho (2004) to model the line profiles 
of the broad and narrow components of Ha and H/3 and use these to measure the total line 
intensity as well as that of the broad components. We begin by fitting the [S II] AA6716, 6731 
profile, when available, fixing the line separations and limiting the line centers to within ±5 
A from the theoretical values. We use up to three Gaussian components to model the line 
profile and allow the scaling between the two [S II] lines to be a free parameter. Using this 
as our narrow-line model, we then fit a narrow Ha+[N II] AA6548, 6583 model fixing the 
profiles of the lines to those we derived from [S II] and the relative strength of A6583 to 
A6548 to the theoretical value of 2.96. We simultaneously fit a broad Ha line allowing up 
to three Gaussian components to accommodate the fit. If [S II] is unavailable, we model the 
narrow-line profile with [O III] A5007 when it has high signal-to- noise. In about a third of 
the objects, no narrow-line profile was available. These are all higher-redshift objects (all 
but one with z > 1) with lower signal-to-noise spectra, especially at H/3. These objects are 
also among the highest luminosity sources (see Figure 12) and have very broad Ha profiles. 
In these sources, we fit up to three Gaussian components to Ha ignoring any (relatively 
weak) narrow-line component. 

To fit a line profile to the broad and narrow components of H/3, we first fit and subtract 
the [O III] AA4959, 5007 doublet, fixing the line separation and allowing the widths and 
intensities to be fit. We perform this subtraction on all our spectra, regardless of whether 
we used the [O III] profile as a narrow-line model, so that even low signal-to-noise emission 
from the [O III] doublet is removed before B.(3 is fit. In the cases where we use a narrow-line 
model, we fit H/3 with narrow- and broad-line components, fixing the widths of the broad-line 
components and the relative positions of the broad and narrow components based on the 
Ha fits. 

We compute the line intensity for the broad-only as well as total line profile, subtracting 
the continuum from the fit. We divide the integrated intensity of Ha by that of H/3 to obtain 
the Balmer decrements. 

To obtain extinction values from our Balmer decrements, we follow Calzetti et al. (1994) 
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In order to make a meaningful comparison between the continuum fit and Balmer decre- 
ment measurements of the reddening, we use the Balmer decrement from the FBQS compos- 
ite spectrum as our intrinsic Ha/H/3 ratio which we measure to be 4.526 for the broad-only 
components and 4.091 for the combined lines, rather than the theoretical value of 2.88 (Os- 
terbrock 1989). Presumably the FBQS composite already includes some reddening; adopting 
the composite Ha/H/3 ratio allows a direct comparison with the E(B — V) values derived 
from the continuum fits to the composite. We also use the same SMC extinction law, k(X), 
that we used in the composite fitting procedure, which gives fc(H/5) — A; (Ha) = 1.259 (using 
the extinction curve parameters from Gordon & Clayton 1998 and the reddening formula 
from Fitzpatrick 1999). 

Figure 11 shows four examples of our Gaussian fitting for computing Balmer decrements. 
The solid line shows the total profile fit. In panels (c) and (d) where narrow-line models 
were used, the Hct+fN II] profile and the narrow H/3 profile are overplotted with a dotted 
line. The broad-line is overplotted with a dashed line. When computing Ha/H/5 we use 
both the total hydrogen line profile as well as only the broad components. We present the 
measurements of E(B — V) based on Balmer decrements from the total line profile listed in 
column (7) of Table 4. 

5.3. Discussion of Reddening 

In figure 12 we plot dereddened absolute i^-band magnitude versus redshift. We color 
the points according to their measured extinctions derived from the continuum fits (from 
unreddened in yellow to heavily obscured in red) and overplot them on top of the FBQS II 
and III quasars which we assume to be largely unabsorbed. The dashed lines indicate the 
detection limit of a flux-limited survey with K < 15.5 for different amounts of extinction. 

The F2M red quasars are among the more luminous objects in this plot especially at 
higher redshifts. Still, even the high-redshift quasars are not luminous enough to allow the 
detection of quasars reddened more than EiB — V) ~ 0.5. At redshifts below z ~ 0.5 we 
find quasars with extinctions of up to E[B — V) ~ 1.5. This trend was noted in White et al. 
(2003b) and in Paper I, and demonstrates an inherent bias in flux-limited surveys against 
finding highly extinguished objects. Even at near-infrared wavelengths, the radiation in the 
rest frame is optical light beyond z ~ 1 and is quite sensitive to extinction. Ongoing quasar 
surveys using the Spitzer Space Telescope should avoid this problem and yield quasar samples 
unbiased by extinction (cf. Lacy et al. 2004; Stern et al. 2004). 

Since almost all the data points in a given spectrum (excluding strong emission lines 
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and regions of strong atmospheric absorption) are used in the composite fitting, this method 
in effect determines the reddening of the continuum. Richards et al. (2003) found that 
intrinsically blue quasars selected from the Sloan Digital Sky Survey may have spectral 
slopes that vary between —0.25 < a u < —0.76. They conclude that these variations are 
intrinsic and not a result of dust extinction. However, our extinction measurements are not 
sensitive to this distinction and may fit an intrinsically "optically steep" quasar with too 
large an extinction value. Therefore, it is useful to compare our derived E(B — V) values 
from continuum fits with E[B — V) derived from Balmer decrements to study possible 
discrepancies between continuum shape and line ratios. 

To determine if the intrinsic slope of a given quasar is responsible for the discrepancy 
between the extinction derived from Balmer decrements and the composite fitting, we varied 
the slope of our composite and recomputed the extinctions. It is well-known that quasar 
spectra, though well-approximated by a power-law over a limited wavelength range, change 
shape in different wavelength regimes and are best described by several broken power-laws 
(Elvis et al. 1994; Zheng et al. 1997; Vanden Berk et al. 2001; Glikman et al. 2006). The 
UV-optical part of the spectrum covering the wavelength range ~ 1300 to ~ 5000 A is 
best described by f u oc z/ -0 - 46 on average (Brotherton et al. 2001; Vanden Berk et al. 2001). 
Therefore we vary the spectral index of our quasar template only over this wavelength range, 
leaving intact the original shape of the rest of the template. To shift the slope to a new 
spectral index, a u , we must multiply the spectrum by a power-law with an index which is 
the difference between the two, Aa u = a u — a u0 . The new template becomes: 



We tested the extreme values found by Richards et al. (2003) (a u = —0.25 and —0.76) 
on all our spectra. From these results we compared the 11 objects for which the Balmer 
decrement fits and continuum fits are discrepant by E{B — V) > 0.2 and whose individual 
fits are very good, since for such cases the discrepancies cannot be attributed to poor data 
and might be explained by an extreme intrinsic slope of the quasar. We use three power-law 
shapes for our template fitting (see Figure 7 in Richards et al. 2003). The default spectrum 
is the mean spectrum derived in Brotherton et al. (2001) which we have used for fitting all 
the red quasar spectra to derive the E(B — V) values listed in Table 4. We also modify 
our template to have the steepest (in wavelength space) spectral shape found by Richards 
et al. (2003) with a u = —0.25. It is clear from this curve that in cases where the Balmer 
decrement measures a higher reddening than the continuum fit, we expect that a template 
with a u = —0.25 will provide better agreement since more extinction is required to redden 




for A < A = 5500 A 



(5) 
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this continuum. When the Balmer decrement underestimates the reddening compared to the 
continuum fitting, it may be because the quasar has a flatter (in wavelength space) intrinsic 
spectrum requiring less extinction to produce the observed spectral shape. 

Four of these objects have E(B — V) > 1 measured from their Balmer decrements. 
Although using the template with a v = —0.25 does yield higher extinctions, the increase 
is typically < 0.1 mag, not enough to match the color-excess derived from the Balmer 
decrements. Similarly, we compared continuum fits using a template with a u = —0.76 to 
four quasars whose Balmer decrements measure little or no reddening (E(B — V) < 0.1). 
Again, though the derived E(B — V) values decreased, it was only by < 0.1 mag. It is possible 
that these objects require even more extreme slopes in their intrinsic spectra; however it is 
unlikely that such a large percentage (8/19) of the objects in our sample with both continuum 
fits and Balmer decrements all have extreme intrinsic continua. 

Figure 13 compares the results of the two techniques for measuring extinction, using the 
same extinction law and with respect to the same intrinsic quasar template. The color-excess 
derived using Balmer decrements of the total hydrogen line profile are plotted with square 
symbols, while the color-excess computed from Balmer decrements of the broad component 
of the lines are plotted with stars. The horizontal error-bars indicate the range of E(B — V) 
values obtained when using templates with the a v = —0.25 and —0.76 spectral slopes. The 
error bars for the E(B — V) from the broad component are much larger than from the total 
line flux, and increase with reddening. A linear fit correlating the two techniques reveals that 
the total line profile (solid line, slope=1.62±0.39, intercept=— 0.35 ±0.23) is closer to a one- 
to-one correlation (plotted with a dotted line for comparison) than the broad component 
(dashed line, slope=2.99 ± 0.57, intercept=— 0.61 ± 0.34). However, even when using the 
total line flux to compute extinctions, this plot shows that although the two techniques 
yield similar reddening measurements on average, the scatter in the fit is ~ 0.5 magnitudes 
and cannot be attributed to a steeper or flatter intrinsic spectrum. Furthermore, Balmer 
decrements can yield drastically higher extinction values by E (B — V) ~ 1 — 2 in noisy 
spectra, since H/5 is a weaker line, especially when reddened, and has large errors in the fit. 

The better correlation between the color-excess derived from Balmer decrements from 
the total line profile, including both the narrow and broad components, with the color- 
excess derived from reddening in the continuum suggests that the dust in these systems is 
not necessarily confined to the "dusty torus" invoked in quasar unification models (Urry & 
Padovani 1995). AGN unification is effective at explaining the Seyfertl to Seyfert2 class 
continuum, where, in this toy-model picture, broad lines seen in the former diminish in the 
latter as the viewing angle approaches 90 degrees from the jet-axis. The narrow-line region 
is farther from the nucleus than the broad line region and is the only contributor of lines 
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in Seyfert2 spectra. Since the Balmer decrements derived from the broad components of 
the hydrogen lines in the F2M spectra correlate poorly with the reddening derived from 
the continuum, it is more likely that the dust in these objects resides outside the nuclear 
AGN, possibly in starforming regions in the host galaxy. Hubble Space Telescope images of 
13 F2M red quasars show and unsusually high fraction of these objects are mergin and/or 
interacting systems. These images suggest that the obscuring dust is broadly distributed 
in these systems (Urrutia et al. 2005). This picture is consistent with the starburst-AGN 
connection suggested by Sanders et al. (1988) and recently modeled by Hopkins et al. (2006). 

In Figure 14 we show the reddening as a function of redshift (for both estimators) using 
all 55 extinction measurements from the continuum reddening (open circles) and the 19 
extinction measurements from the Balmer decrements determined from the total line fluxes 
(filled circles). The error bars on the open circles indicate the range of E(B — V) values 
obtained from fitting the templates with extreme steep and flat spectral slopes. We connect 
the points from the same object with a line. The discrepancies between the two techniques 
do not appear to depend on redshift. However, both techniques do show the bias toward 
redder objects at lower redshifts noted above. 

5.4. Are There Variations in the Dust Law? 

In addition to the variations in the quasar template, variations in the dust law could 
also affect our color-excess measurements and the quality of our fits. Figure 3 shows that the 
F2M red quasar sample includes a highly non-uniform set of spectra. In addition to intrinsic 
variations in the spectral shapes of these quasars, our wavelength sampling and coverage as 
well as signal-to-noise are also highly diverse. In cases where the continuum fits have not 
been improved by changing the slope of the quasar template, it is possible that the dust law 
itself is different enough to affect the fit. 

The extinction law for the Milky Way Galaxy has been studied extensively using UV- 
though-infrared observations of stars along many lines of sight. The extinction as a function 
of wavelength is typically parameterized by the ratio of total to selective extinction in the 
V-band, R v = A(V)/E(B - V) (Cardelli et al. 1989). Galactic values range from 2.2 to 
5.5 with a mean value of ~ 3.1. Our extinction law is that derived by Gordon & Clayton 
(1998) for the Small Magellanic Cloud and is shown in Figure 15. The solid line represents 
Rv = 3.1, while the more extreme values are overp lotted. The only difference between the 
dust laws occurs in the near-infrared. Dust laws for dust composed of larger grains are 
characterized by larger values of Ry and have a weaker effect in the near-infrared causing 
less extinction. Shortward of ~ 7000A, the dust laws behave the same, and therefore we 
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do not expect differences in the dust law to affect the extinctions that we derive from the 
Balmer decrements. 

We computed E(B — V) for all 55 of our quasar spectra using three extinction curves, 
with R v = 2.2, 3.1, and 5.5, and compared the results. The quality of the fits were mostly 
indistinguishable. In a few cases where the combined optical plus near infrared spectra were 
well-fit by the original Ry = 3.1 dust law, the extreme values produced slightly poorer fits. 
The fits to the Ry = 2.2 dust law yielded higher E(B — V) values while the Ry = 5.5 law 
yielded lower E(B — V) values. Figure 16 plots the color excess as a function of redshift 
for Ry = 3.1 in filled circles. The error bars represent the extreme extinctions derived 
from the Ry = 2.2 and 5.5 dust laws. The results varied most in spectra for lower redshift 
objects where we only had a near infrared spectrum. In these cases the varying sensitivities 
of the dust laws produced strongly differing results, since the Ry = 2.2 dust law is more 
sensitive in the infrared and therefore a higher E(B — V) value is computed. When only 
an optical spectrum was available, or at higher redshifts where rest frame optical light is 
in the observed near-infrared spectrum, the derived extinctions were the same regardless of 
the dust law. Combined spectra were mildly affected, and only at low redshifts. Therefore 
the apparent trend with redshift is more a result of the heterogeneous nature of our spectra, 
than of any physical effect of dust law with redshift. To better determine the nature of the 
dust in our spectra, longer wavelength data are needed (e.g., Spitzer IRS spectra). 



6. Radio Spectral Indices 

We have so far explored the red colors of our F2M quasars by assuming that they are a 
result of dust reddening. It is possible that there is also a contribution to the i^-band flux 
from nonthermal synchrotron emission from jets. Serjeant & Rawlings (1996) argued that 
enhanced synchrotron emission from relativistic jets lying in the line of sight contributes to 
the near-infrared flux and is responsible for both the flat radio spectrum and the red optical 
to near-infrared colors in the red quasars found by Webster et al. (1995). Whiting et al. 
(2001) looked for a synchrotron signature in the spectral energy distributions of the Parkes 
red quasars and found evidence for a red synchrotron component in the near-infrared. 

The Parkes quasars were drawn from the Parkes Half-Jansky Flat-spectrum Sample 
(PHFS; Drinkwater et al. 1997) which included only sources brighter than 0.5 Jy at a radio 
frequency of 2.7 GHz. It is possible for these synchrotron spectra to extend toward higher 
frequencies and contribute enough near-infrared flux to account for the large R — K colors. 
On the other hand, the FIRST survey is sensitive down to 1 mJy and 43/56 of the F2M 
quasars are fainter than 30 mJy at 20 cm (1.4 GHz). We can put limits on the K-band 
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light expected from synchrotron radiation by extrapolating the radio spectrum into the near 
infrared. Figure 17 shows the magnitude of a series of 1 mJy FIRST sources (y-axis on the 
right) with flat to inverted spectral indices (a u ranging from —0.5 to 0.2) at 2.2 /xm (K band, 
y-axis on the left). This computation makes the simple (and probably naive) assumption 
that there are no breaks in the radio spectrum. Although this is likely to be unrealistic, it 
gives us the "worst case scenario" of how much synchrotron flux could be contaminating the 
K band part of the spectrum. Under these assumptions, a radio source with a flux density 
of 1 mJy at 20cm with a radio spectral index of > 0.0 could contribute < 14.6 magnitudes in 
K band, which would be detectable in 2MASS. A 10 mJy source at 20 cm would contribute 
14.6 magnitudes in K band even with a spectral index of a u = —0.2 and 12 magnitudes with 
a„ = 0.0. 

Since the faintness of our radio sources does not rule out the possibility that a flat- 
spectrum source can contaminate our K band light and redden otherwise blue, unobscured 
quasars, it is necessary to check whether the F2M quasars (selected to be red) have a higher 
fraction of flat-spectrum sources (accounting for their color) compared with optically-selected 
radio-detected quasars, such as FBQS. 

The quasars in our study reside in a portion of color-color space largely unexplored 
by optical surveys. Of the more than 14,000 quasars in the Veron-Cetty & Veron (2001) 
catalog, there are only 21 with 2MASS matches that fall outside our survey area, but that 
reside in our prescribed region of color-color space (R — K > 4 and J — K > 1.7). Nineteen 
of the 21 are radio-selected objects. All of these have flat spectral indices (a > —0.5) when 
radio fluxes at two or more frequencies are available. 

These findings are suggestive but are not conclusive. Furthermore, they do not probe 
sources of intermediate radio loudness, with 20 cm flux densities in the 1-30 mJy range, such 
as the majority of the objects in the F2M survey. Since the fraction of quasars that are 
radio sources rises steeply with a decreasing flux density threshold (White et al. 2000), it is 
important to see if a large number of missing (red) quasars continues to appear as we go to 
fainter fluxes. 

Webster et al. (1995) and others have shown that quasars selected by their flat radio 
spectra have redder colors than optically selected quasars. In order to conclude that the two 
phenomena are linked, optically selected red quasars ought to have flat radio spectra as well. 
To test this, we obtained simultaneous 3.6 cm and 20 cm flux density measurements for 44 
of our 57 F2M quasars using the VLA in the CnD configuration in order to determine their 
radio spectral indices and to test whether the F2M sample exhibits enhanced synchrotron 
emission. Data were reduced in the standard manner using the AIPS package. 
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Flux densities were derived by fitting two-dimensional elliptical Gaussians to the images; 
for the weakest sources, the Gaussian peaks were fixed at the optical positions and the source 
sizes were fixed to the size of the synthesized beam. Peak and integrated flux densities were 
derived; since we are interested in the contribution any core component makes to the IR 
emission from the nucleus, we use the peak flux densities to calculate the spectral indices 
(all but 2 of the 42 objects detected at 3.6 cm are point-like). 

The results are presented in Table 6. We list the quasar's name in column (1), followed 
by the FIRST peak and integrated flux densities in columns (2) and (3), respectively. Our 
newly measured 20 cm peak and integrated flux densities in the VLA CnD configuration 
are listed in columns (4) and (5), respectively. In column (6) we list the NVSS 20 cm flux 
density in the VLA D configuration as measured by Condon et al. (1998). Columns (7) and 
(8) list our newly measured 3.6 cm peak and integrated flux densities, again in the CnD 
configuration. Column (9) lists the radio spectral-index, cci.4 ghz/8.3 (GHz), computed from 
the peak flux densities in columns (4) and (7). 

We compare these measurements with contemporaneous 20 cm and 3.6 cm VLA mea- 
surements for 214 FBQS quasars (Lacy et al. 2001). We divide the samples into a radio-bright 
and radio-faint sample by their 20 cm flux density to avoid any biases in the FBQS sub- 
sample which is biased toward brighter radio sources and thus has a higher proportion of 
flat-spectrum objects. Figure 18 presents spectral index histograms for the F2M and FBQS 
quasars, normalized by the number of quasars in each sample for ease of comparison. The 
bright samples (left panel) have indistinguishable distributions. There is a large fraction 
of flat-spectrum objects (cci.4ghz/8.3GHz > —0.5) in this subset and enhanced synchrotron 
emission should be considered as a possible source for the red colors in these quasars. The 
faint sample (right panel) shows an excess of steep spectrum quasars in the F2M sample. 
This may be due to the slightly lower median 20 cm flux density (2.68 mJy for F2M versus 
3.1 mJy for FBQS). In these objects dust is the likely cause for the reddening. 

In particular, in only eight of the 44 objects observed would a smooth extrapolation 
from 3.6 cm to 2/xm significantly boost the K-band magnitude. Three of these eight are 
blue quasars (see §6), and so are not considered part of the red sample. The other five 

(representing less than 10% of our red quasar sample) have E{B — V) values of ~ 0.2 0.7, 

intermediate values for this population. Thus, we conclude that synchrotron contamination 
does not have a significant impact on any of our conclusions. 
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7. The Intrinsic Surface Density of Quasars on the Sky 

Having measured the extinctions for our red quasars, we have corrected their observed 
.fT-band magnitudes to their unobscured values. The amount of i^-band absorption ranges 
from Ak ~ — 1.2 and once corrected, the intrinsic magnitude-limit of our survey becomes 
K < 15.0. Figure 19 shows the density of quasars as a function of dereddened ^-magnitude. 
The number of red quasars begins to decline at K > 14.5. This is a result of the 2MASS 
limit of K ~ 15.5 corrected for a median (A K ) ~ 0.5. We ignore the 36 candidates without 
spectra, among which are the 19 objects that are undetected in the 2MASS J-band with 
J — K < 1.7 (half of which have J - K < 1.4). We also omit F2MJ021542.0-022257, the 
z = 1.187 quasar from Drinkwater et al. (1997), since we do not have a spectrum from which 
to measure E(B — V) and determine its intrinsic i^-band magnitude. These objects were 
included in Figure 7. 

We can assume our quasar detection efficiency of 50% applies to the 17 objects with 
J-band detections and expect ~ 8 additional red quasars whose observed magnitude range 
lies between 14.6 < K < 15.5 (the fainter half of the /A'-magnitude distribution for the 
whole sample, see Figure 4). Since it is possibile that many of the objects lacking a J-band 
detection actually have J — K < 1.7 and are therefore likely to be stars, the probability that 
these objects are quasars is likely less than 50%. They fall into the same 14.6 < K < 15.5 
magnitude range as the former group. 

For all the unclassified candidates, we estimate that there are fewer than 18 quasars. 
Assuming the absorption for these objects is the median A K of the rest of the quasar sample 
puts them in the dereddened range of 14.1 < K < 15. These objects may increase the 
density of objects with dereddened K > 14.0, where our sample appears to be significantly 
incomplete. Of course, their extinctions may be higher than the median, adding counts to 
the density measurements at K < 14.0. 

Below K = 14.0 we are able to make a comparison between the dereddened F2M quasars 
and FBQS II and III (with poorer statistics). The density of F2M quasars with dereddened 
K < 14.0 is 0.011 ± 0.002 deg" 2 . At the same limit, FBQS II has a density of 0.028 ± 0.003 
deg -2 (summing the first six rows in column (1) of Table 5), implying a missed red quasar 
fraction of 28 ±6%, using FBQS II. The FBQS III survey covers a smaller area and therefore 
has fewer objects brighter than K = 14. The density of FBQS III objects at this limit is 
0.008 ± 0.004 deg -2 , implying a missed red quasar fraction of 56 ± 16%. 

To a dereddened limit of K < 14.5 the spatial density of F2M quasars is 0.016 ± 0.002 
deg -2 which implies that red quasars make up 22 ± 4% of the total quasar population, using 
the FBQS II sample, and 30 ± 5% using FBQS III. These values are lower limits, since the 
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17 unobserved quasars may contribute counts to the F2M sample at these magnitudes; in 
addition, quasars with (A K ) > 1.2 would fall below the 2MASS threshold. 

Again, these are likely lower limits, based on the arguments put forth in Section 4. In 
addition, there may be redshift-dependent magnitude biases that exist for reddened objects 
(as we discussed in Section 5.3) which may impact our estimates of the fraction of red 
quasars. Our plot of dereddened absolute i^-magnitude versus redshift (Figure 12) shows 
regions where the red quasars appear to make up more than 25 — 50%. For example, the 
Mk < —30, 0.5 < z < 1 region appears to have about 7 red quasars and 6 unreddened 
quasars, before making any correction for red objects that can't be detected. On the other 
hand, the number of red quasars at the fainter, lower redshift end, where the selection effects 
against finding red quasars are weaker, appears significantly smaller than 50%. To untangle 
the combined effects of luminosity, reddening and redshift on a flux-limited survey, even in 
the K band, requires a larger sample reaching fainter flux-limits. 



8. Blue Quasars in our Red Quasar Survey 

Though our survey is highly efficient at finding red quasars by selecting for red colors 
in two sets of filters, we find a handful of blue quasars contaminating our survey. In Paper I 
we found one blue quasar, F2M J090651. 5+495235, at a redshift of 1.635. We recover this 
source, along with four others in this survey. Why do we find blue quasars in a survey 
targeted for red objects? 

We define "blue" quasars as those with E(B — V) < 0.1 (in column 8 of Table 4). In 
Tables 7 and 8 we present data we have collected from the literature on these objects. We 
present optical photometry from three different epochs, GSC II, APM and SDSS to look 
for evidence of variability to explain why blue quasars end up in our survey; e.g., objects 
observed at a bright phase with 2MASS but a faint phase in the POSS-II observations. We 
also present the FIRST and NVSS radio flux densities to search for variability in the radio 
which may further suggest variability in the optical. Four of the five (80%) blue quasars have 
20 cm flux densities above 100 mJy, compared with 6 of the 52 (12%) red quasars in the whole 
F2M sample 5 ; since radio-loud quasars are more highly variable than radio-intermediate and 
radio-quiet populations, this discrepancy is consistent with variability being the explanation 



5 We compare radio brightness, rather than radio-loudness (defined in Table 8) because the large extinction 
in the optical will produce an artificially large radio-loudness. To compute the true value of radio-to-optical 
ratio we must de-redden the B magnitudes of our quasars by Ab, which, as we have shown in §5, can be 
highly uncertain. 
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for at least some of the blue quasars in our sample. 

Tables 7 and 8 summarize the optical and radio properties of these objects, respectively. 
Columns (1) and (2) of Table 7 list the J2000.0 Right Ascension and Declination, followed 
by the redshift in column (3). The optical photometry is presented in the following nine 
columns. Columns (4), (5), and (6) present the GSC-II catalog's photographic F (red) and 
J (blue) magnitudes, and J — F color, respectively. Columns (7), (8), and (9) list the APM 
photographic E (red) and O (blue) magnitudes, and O — E color, respectively. Columns 
(10), (11), and (12) list the SDSS r and g magnitudes, and g — r color, respectively. Column 
(13) lists the 2MASS i^-band magnitude and column (14) lists the color excess, E(B — V), 
derived from Table 4. Table 6 presents the J2000.0 coordinates of the five blue quasars 
in columns (1) and (2) followed by the FIRST peak and integrated radio flux densities in 
columns (3) and (4). We list the NVSS 20 cm flux density in column (5) and the 6 cm 
flux density, mainly derived from Becker et al. (1991), in column (6). We present the radio 
spectral index, a20cm/6cm, using the 20 cm flux density from the FIRST (integrated) and the 
NVSS flux density in columns (7) and (8), respectively. Column (9) lists the radio- loudness 
parameters. Column (10) indicates whether the object is a well-known radio source, based 
on references listed in the NASA/IPAC Extragalactic Database (NED) and column (11) 
lists any comments on the source based on observations at other wavelengths or from the 
literature. 

All of these objects have flat or inverted radio spectra (a > —0.5, S u oc u a ) and have 
been studied at multiple wavelengths. F2M115931.8+291443 and F2M164258.8+394837 are 
well-known Blazars with well-studied variability properties from 7-rays to radio wavelengths. 
F2M115931.8+291443 (4C +29.45) has shown optical variability of AV > 5 magnitudes 
(Wills et al. 1983; Ghosh et al. 2000). F2M164258.8+394837 (3C 345) underwent an outburst 
in 1991 of AB ~ 2.5 magnitudes (Webb et al. 1994; Raiteri et al. 1998). A similar change of 
~ 2 magnitudes is seen between the APM and GSC-II photometric data presented in Table 
7. 

Figure 20 shows our three other blue quasars overplotted by their SDSS g, r, i and 
2MASS J, H, K s photometry. The objects in the top two panels show a strong discontinuity 
between the optical photometric points from SDSS and the 2MASS near-infrared photometry. 
These objects also have inverted radio spectral indices consistent with radio and possibly 
optical variability. The spectrum in the bottom panel shows a discontinuity between the 
optical and near-infrared photometry, but also indicates a flattening in the optical spectral 
shape compared with the spectrophotometric points. 

Giommi & Colafrancesco (2004) reveal that the spatial density of Blazars (also referred 
to as Flat-Spectrum Radio Quasars; FSRQ's) brighter than 0.1 Jy at 5 GHz is 0.18 deg -2 . 
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The spatial density of featureless BL Lac objects brighter than 0.1 Jy at 5 GHz is 0.025 
cleg" 2 . 

Our survey is 77% complete over 2716 deg 2 suggesting 376 (2716 deg 2 x 0.77 x 0.18 
deg~ 2 ) expected Blazars and 52 (2716 deg 2 x 0.77 x 0.025 deg" 2 ) BL Lacs. All five sources 
with computed spectral indices in Table 8 are have S^ghz ^ 0.1 Jy. Even if we assume that 
all five of these quasars are Blazars (two are confirmed as such in the literature and the 
other three have strong circumstantial evidence that they belong to this class) and add the 
five BL Lacs from the main sample, we are finding less than 10% of the numbers expected, 
demonstrating that our survey for red objects strongly selects against these sources. 

9. Red Quasars in Radio Surveys 

In Section 6 we showed that the red colors of the F2M quasars are unlikely to be caused 
by synchrotron radiation and are most probably a result of dust extinction. Although it is 
hard to conceive of a causal connection between the radio properties of quasars and their 
extinction, there are far more radio-emitting red quasars known than radio-quiet red quasars. 

The best explanation for this is that reddened quasars elude optically-selected samples 
both because they fall out of flux-limited samples and because their colors begin to resemble 
those of low mass stars, a huge source of contamination. Samples selected in the radio are 
insensitive to this bias. Although only 10 — 20% of quasars are radio-emitters, stars are 
effectively radio-silent. Radio selected quasar samples are therefore more effective at finding 
red quasars. And quasar composites made from radio-loud and radio-intermediate quasars 
tend to have redder colors than radio-quiet quasars for this reason (cf. Brotherton et al. 
2001). 

There are several examples of obscuration in very radio-bright samples, such the Third 
Cambridge (3C) catalog of radio sources (whose members are brighter than 10 Jy). Hill 
et al. (1996) looked for broad Paa in low-redshift narrow-line radio galaxies (NLRGs) from 
the 3C radio catalog. They reclassify three of the NLRGs as broad-line radio galaxies after 
detecting broad Paa in their near-infrared spectra, as was the case with some of the F2M 
red quasars. Rawlings et al. (1995) found a broad Ha line redshifted into the near-infrared 
in 3C22, whose observed optical spectrum shows only narrow lines. Similar to many of the 
F2M red quasars, 3C68.1 shows broad lines in the optical atop a reddened, highly polarized 
continuum (Brotherton et al. 1998). Without their strong radio emission, none of these 
objects would have been selected as quasars in optical surveys. 

Hill et al. (1996) were able to explain the extinction in their radio galaxies in terms of 



-28- 



orientation. They found that the broad lines were more reddened than the narrow lines in 
their spectra, placing the absorbing dust between the broad- and narrow-line regions. We 
were not able to confirm this finding for the F2M quasars. In fact, as we showed in Section 
5.3 and Figure 13, the reddening that we measured from the continuum correlated more 
strongly with the reddening derived from Balmer decrements that were measured from the 
total (broad plus narrow) line profiles, rather than the broad-only components. This suggests 
that the dust in the F2M quasars is probably dispersed throughout the host galaxy. The 
obscuration seen in our sources is likely not an orientation effect, but rather an evolutionary 
one. 

Although the F2M red quasars are radio-emitters, by definition, they are not neces- 
sarily radio-loud. It is misleading to apply the traditional definition of radio-loudness, 
R = /(5GHz)//(B) (Stocke et al. 1992), to the F2M quasars without correcting for ex- 
tinction, as most (70%) would fall into the radio-loud regime. However, to compute the 
unreddened 5-band flux-density we must correct for the extinction, which, as we have ar- 
gued in Section 5.1, is highly uncertain. Paper I computed R for the optically-faintest subset 
of the F2M sample and found that two-thirds of the quasars fell into the radio-intermediate 
regime (3 < R < 100), similar to the FBQS quasars. 



10. Summary and Conclusions 

We have presented the FIRST-2MASS red quasar survey based on sources matched 
between the FIRST and 2MASS catalogs with optical-to-near-infrared colors obeying J— K > 
1.7 and R — K > 4. Our catalog consists of 156 quasar candidates of which 120 have been 
spectroscopically identified. We have found 57 quasars, ~ 50 of which are true red quasars, 
most likely obscured by dust. 

We have applied two techniques for determining the extinction in these sources and have 
compared the results. By applying a Small Magellanic Could extinction law to a quasar 
template, we fitted the spectral shape of our reddened quasars and derived E(B — V). 
In objects where Ha and Hj3 appear in the same (optical or near-infrared) spectrum we 
compute extinctions based on the Balmer decrement for both the broad-only component of 
the lines as well as the total Ha and H/3 line profiles. We find better agreement between 
the two techniques when comparing to the Balmer decrements derived from the total line 
profile, though there is a large scatter (~ 0.5 magnitudes) between the extinction values 
derived from the two techniques. We investigate the effect of changing the optical slope of 
the quasar template and find that this has a small impact on the values derived. We also 
explore varying the grain size in our dust law by changing the Ry parameter. The results 
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from these tests were inconclusive due to the highly non-uniform nature of our spectra. Since 
the SMC dust law varies only at long wavelengths (see Figure 16), SpitzerlRS spectra would 
cover the ideal wavelength regime to determine the dust law for these quasars. 

We note a bias in our sample based on the flux-limited nature of our survey (and lower 
spectroscopic completeness at fainter magnitudes) which prevents our detection of highly 
reddened quasars at high redshift. Our sensitivity to reddened quasars as a function of 
redshift declines and we can only detect the most luminous sources with moderate extinction 
at z > 1.5. Deeper K-b&nd surveys such as the UKIDSS survey (sensitive to K — 18.4) will 
provide an opportunity to detect highly obscured quasars at higher redshifts. We can then 
begin to look for dependence of extinction on redshift and intrinsic luminosity. 

We have estimated the percentage of missed red quasars in flux-limited optical surveys 
that would be recovered in flux-limited near-infrared surveys. We find that optically selected 
flux-limited quasar surveys (such as the FBQS II and III) miss ~ 12% of quasars as a 
consequence of extinction pushing the observed magnitudes below the surveys' optical flux 
limit. We also estimate the true fraction of red quasars after correcting for extinction effects. 
We find that for quasars brighter than K = 14.0 red quasars make up between ~ 25% 
and 60% of the total quasar population (depending on the optically selected comparison 
sample). Below K = 14.5 we find that red quasars make up > 20% to 30% of the total 
quasar population; this represents a lower limit due to our survey's incompleteness at fainter 
K-band magnitudes. 

Finally, we looked at the handful of blue quasars that we find in this survey, concluding 
that most are selected as a result of their optical variability. We present radio spectral indices 
for a subset of our red quasars based on contemporaneous 20 cm and 3.6 cm measurements 
with the VLA. We conclude that synchrotron radiation may only contribute to the quasars' 
red colors in the most radio-bright sources (S20 cm > 10 mJy), but the distribution of radio- 
spectral indices of our sample suggests that reddening as a result of dust obscuration cannot 
be explained by orientation models for the F2M quasars with flat radio spectra. 
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07 


G104 
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1.53 


1.45 


20.81 
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16.610 


15.475 


14.836 


1.77 


4.03 


SbGal 
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09 21 22.13 


+30 29 09.9 


1.25 


2.15 


21.59 


20.38 


17.582 


16.192 


14.993 


2.59 


5.39 


Galaxy 


0.423 


O10 




09 21 45.69 


+19 18 12.6 


4.40 


4.72 


21.32 


20.12 


16.749 


15.429 


14.548 


2.20 


5.57 


QSO 
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09 IR8 
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4.73 
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15.465 
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16.933 


15.047 
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Galaxy 
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8.76 
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15.270 


2.06 


>5.53 


QSO 
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19.56 
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NLAGN 
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G104 
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66.85 


20.70 
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17.811 
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1.00 
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15.124 
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SbGal 
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2.09 
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15.121 
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22.37 


21.19 


20.42 
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4.74 
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15.888 


14.925 


1.97 


4.05 










10 43 35.72 


+19 28 21.5 


1.16 


1.44 


21.94 


19.62 


17.630 


16.502 


15.552 


2.08 


4.07 


Galaxy 


0.570 


Oil 




10 43 55.37 
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1.74 


1.78 


21.32 


19.65 


17.713 


16.590 


15.488 


2.22 


4.16 


Galaxy 


0.399 
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02 
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10 


31 


.6 
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15 
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2 
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10 49 


18 
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14 


58 


.1 
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15 


.362 


2 


49 
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07 
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+28 


12 


47 


.0 


211.47 
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16. 


.556 
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544 


14 
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1 
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11 13 


54 
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14 


38 


.9 


2.99 


4.22 


20.87 


19.00 


16. 
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15. 
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13 
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2. 


.55 
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25 


.05 
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55 


30 


.4 


5.51 


6.57 


21.51 


19.24 


17 
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15 


.717 


15 


.238 


2 


00 


11 17 


33 


.83 


-02 


36 


00 


.3 


608.04 
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20.27 
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15, 


.952 


14. 


.980 


14 
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1. 


90 
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11 
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33 


41 


.9 


1.30 


1.81 


>22.50 
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14 
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2 
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19 
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55 


31 
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1 
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56 
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3 
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58 
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58 


34 


.2 
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13 
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1 
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11 51 


24 
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59 


57 


.4 


3.52 


3.60 


21.92 
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17. 


105 


16. 
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15 
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2. 


.00 


11 53 


14 


.48 


+29 


05 


23 


.8 


4.66 


4.81 


21.68 


19.25 


16. 


.903 


16 


.248 


15 


.119 


1 


.78 


11 59 


31 


.84 


+29 


14 


43 


.9 
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17.67 
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13 
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12. 


.304 


11 


.477 


1 
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12 01 


11 


.14 


-03 


32 


19 


.7 


55.27 


56.56 


>22.50 


19.40 


16. 


.680 


15 


.849 


14 


.921 


1 


.76 
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55 


.33 
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15 


18 


.8 


43.01 
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19.34 


17 


.276 


16 
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15 
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2 
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12 06 


18 
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57 


19 


.7 
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5.59 
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16 
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15 
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2. 
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21 
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07 


17 


.0 
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1.25 
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15 
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13 
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2. 


20 
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.2 
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02 


12 27 


03 
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56 
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4.27 
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110.98 
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15 


.582 


14 


.264 


3 


.15 


13 57 


32 


.14 


+32 


45 


29 


.8 


1.28 


3.48 


21.73 


19.36 


17 


.174 


16 


688 


15 


138 


2 


.04 


13 59 


41 


.18 


+31 


57 


40 


.5 


1.33 


1.28 


>22.50 


>20.80 


16. 


927 


15. 


835 


14 


.788 


2. 


,14 


14 08 


11 


.61 


+32 


43 


50 


.2 


1.09 


2.26 


>22.50 


19.59 


17. 


.152 


16 


.168 


15 


.353 


1 


80 


14 12 


14 


.30 


+29 


12 


08 


.3 


14.34 


14.22 


20.87 


19.42 


17. 


.287 


15. 


.499 


14 


.890 


2. 


,40 


14 15 


22 


.83 


+33 


33 


06 


.5 


5.21 


6.21 


20.58 


18.92 


16 


.716 


15 


.372 


14. 


.304 


2 


,41 



4.47 
4.73 

>5.44 
5.39 
5.33 
4.00 
4.70 
5.29 
4.01 
6.06 
4.64 
>5.70 
4.13 
6.11 
4.48 
4.15 
4.03 
4.14 
4.20 
4.16 
4.39 
4.00 
5.92 
4.62 
>5.95 
>6.02 
>6.54 
4.22 
>6.01 
4.24 
4.53 
4.62 



QSO 


1.440 




X94 


SbGal 


0.561 


07 




NLemit 


0.715 


07 




Galaxy 


0.677 


015 


G104 


QSO 


0.842 






QSO 


0.681 


05 IR1 




Galaxy 


0.433 






QSO 


0.463 


04 




QSO 


0.686 


05 IR8 


G104 


NLAGN 


0.394 


06 




QSO 


1.088 


Oil 




QSO 


0.780 


05 IR1 


G104 


QSO 


0.730 




■ 

Gr96 


Galaxy 


0.755 


013 




NLAGN 


0.566 


05 


G104 


Galaxy 


0.166 


08 




QSO 


0.363 


IR5 








021 




QSO 


0.768 


020 IR6 




BLLac 


0.993 






QSO 


0.315 


IR5 




QSO 


0.275 


05 




NLemit 


0.210 


IR3 




QSO 


1.720 


Oil IR3 


G104 


QSO 


1.770 


Oil IR3 




QSO 


1.311 


09 IR8 


G104 


QSO 


1.724 


IR1 


G104 


QSO 


0.338 


012 




QSO 


0.416 


IR5 





CO 

GO 



Table 1 — Continued 



R.A. 


Dec 






B 


R 


J 


H 


K s 


J-K 


R — K 


Type 


Redshift 


Observations 


Ref. 


(J2000) 


(J2000) 


(mjy) 


(mjy) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 










(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


14 22 46.15 


+24 04 11.4 


1.87 


1.35 


>22.50 


>20.80 


17.525 


16.168 


15.325 


2.20 


>5.47 


Mstar 




Ol 


G104 


14 25 49.04 


+14 24 56.7 


509.19 


523.79 


19.28 


19.19 


16.636 


15.941 


14.843 


1.79 


4.44 










14 27 44.34 


+37 23 37.4 


1.71 


1.52 


21.50 


19.59 


17.110 


15.780 


15.071 


2.04 


4.52 


QSO 


2.168 


Oil 




14 46 59.92 


-00 50 14.9 


30.10 


35.69 


21.52 


18.94 


16.745 


16.080 


14.801 


1.94 


4.14 






IR6 




14 48 12.15 


+30 56 14.5 


1.02 


0.94 


22.02 


19.39 


16.672 


16.197 


14.960 


1.71 


4.43 


Galaxy 


0.271 


Oil IR1 


G104 


14 54 06.68 


+19 50 28.0 


1.34 


0.86 


20.74 


18.89 


17.079 


15.915 


14.710 


2.37 


4.18 


Seyl.5 


0.260 




S02 


14 58 43.41 


+35 42 57.5 


688.50 


693.13 


20.63 


18.96 


16.795 


15.888 


14.951 


1.84 


4.01 






019 




14 58 44.82 


+37 20 21.6 


266.39 


269.98 


18.63 


>18.50 


16.128 


15.141 


14.269 


1.86 


>4.23 


Galaxy? 


0.333 




V96 


15 31 50.48 


+24 23 17.6 


2.24 


2.00 


20.52 


19.17 


16.667 


15.783 


14.693 


1.97 


4.48 


QSO 


2.287 


IR6 




15 32 33.19 


+24 15 26.8 


7.43 


9.37 


20.57 


20.09 


16.831 


15.658 


14.956 


1.87 


5.13 


QSO 


0.562 






15 45 16.90 


+38 15 39.6 


1.86 


1.33 


19.96 


18.20 


16.612 


15.314 


14.069 


2.54 


4.13 


NLemit 


0.290 


IR5 




15 49 38.80 


+12 45 07.9 


1.80 


1.79 


18.92 


17.61 


15.824 


14.565 


13.524 


2.30 


4.09 


QSO 


2.373 


IR5 




15 55 23.72 


+34 19 39.1 


1.49 


1.27 


20.84 


19.21 


17.136 


16.151 


15.136 


2.00 


4.07 








1 


16 00 34.56 


+35 22 27.0 


3.17 


3.30 


20.52 


18.59 


16.412 


15.400 


14.170 


2.24 


4.42 


QSO 


0.707 


IR4 


CO 

> 

| 


16 12 53.62 


+33 20 10.3 


2.91 


2.32 


21.08 


19.39 


17.417 


16.130 


15.347 


2.07 


4.04 








16 14 54.10 


+35 14 13.2 


1.87 


1.43 


22.12 


19.32 


17.167 


16.136 


15.306 


1.86 


4.01 










16 18 09.72 


+35 02 08.5 


13.73 


16.22 


>19.50 


18.25 


16.785 


15.408 


14.054 


2.73 


4.20 


QSO 


0.447 


IR5 




16 42 58.80 


+39 48 37.2 


6050.06 


6598.61 


17.42 


16.68 


14.260 


13.289 


12.360 


1.90 


4.32 


QSO 


0.593 


020 IR4 




16 54 07.69 


+34 35 27.4 


1.85 


1.72 


19.97 


18.27 


16.548 


15.386 


14.247 


2.30 


4.02 


NLAGN 


0.215 






16 56 47.11 


+38 21 36.7 


4.12 


3.89 


>22.50 


>20.80 


17.299 


17.718 


15.087 


2.21 


>5.71 


QSO 


0.732 


02 


G104 


17 04 25.31 


+30 37 08.7 


1.33 


1.22 


20.52 


>20.80 


17.311 


16.263 


14.917 


2.39 


5.88 


QSO 


0.311 


IR4 




17 08 02.53 


+22 27 25.7 


1.31 


1.74 


22.32 


19.35 


17.036 


15.959 


14.575 


2.46 


4.77 


QSO 


0.377 


03 IR6 




17 10 58.40 


+31 09 58.1 


6.09 


6.16 


16.32 


>20.80 


16.742 


15.818 


14.853 


1.89 


>5.95 






IR7 




17 24 05.43 


+40 04 35.9 


380.70 


456.24 


21.90 


19.75 


16.130 


15.001 


14.074 


2.06 


5.68 


Galaxy? 


1.49 


IR3 IR5 


V96 


17 24 28.03 


+33 31 19.9 


1.58 


1.35 


20.88 


18.48 


16.927 


15.734 


14.461 


2.47 


4.02 


QSO 


0.370 






17 33 53.18 


+41 24 55.5 


1.63 


0.96 


21.03 


19.04 


17.226 


15.632 


15.003 


2.22 


4.04 










22 16 33.73 


-00 54 51.2 


1.31 


0.80 


19.92 


17.91 


15.946 


14.942 


13.777 


2.17 


4.13 


QSO 


0.200 


018 IR2 




22 22 52.78 


-02 02 57.4 


2.43 


2.66 


21.37 


20.04 


17.436 


16.247 


15.141 


2.30 


4.90 


QSO 


2.252 


010 IR2 




22 24 38.35 


-00 07 50.9 


1.41 


1.29 


>22.50 


19.78 


16.810 


15.703 


15.020 


1.79 


4.76 


Galaxy 


0.452 


O10 


G104 


22 47 30.20 


+00 00 06.7 


183.71 


190.56 


17.99 


18.17 


15.649 


14.820 


13.860 


1.79 


4.31 


BLLac 




019 IR2 


B01 


23 11 04.66 


-02 14 58.2 


4.41 


7.05 


22.05 


19.51 


17.281 


15.588 


15.055 


2.23 


4.45 


Galaxy 


0.499 


07 




23 25 49.88 


-10 52 43.4 


1.46 


1.10 


20.71 


19.70 


17.382 


16.311 


15.238 


2.14 


4.46 


QSO 


0.564 


O10 IR2 




23 28 21.74 


-11 07 03.0 


6.13 


5.82 


19.25 


17.96 


16.051 


15.028 


13.912 


2.14 


4.05 


NLAGN 


0.232 


014 





Table 1 — Continued 



R.A. 


Dec 






B 


R 


J 


H 


K a 


J - K 


R-K 


Type 


Redshift 


Observations 


Ref. 


(J2000) 


(J2000) 


(mjy) 


(mjy) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 










(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


23 39 03.83 


-09 12 21.0 


4.34 


4.06 


20.95 


18.63 


15.805 


15.052 


14.088 


1.72 


4.54 


QSO 


0.660 






23 45 54.91 


-10 03 29.1 


2.39 


2.09 


21.18 


18.48 


16.893 


15.754 


14.187 


2.71 


4.29 


Galaxy 


0.263 


06 IR2 




23 50 40.11 


-10 43 59.2 


1.35 


0.96 


21.27 


19.24 


17.096 


15.918 


15.043 


2.05 


4.20 


Galaxy 


0.253 


IR4 





Note. — Optical spectra: (Ol) ESI 2000 July; (02) ESI 2000 August; (03) ESI 2001 March; (04) ESI 2001 May; (05) ESI 2001 June; (06) ESI 2002 January; (07) 
ESI 2002 February; (08) ESI 2002 March; (09) ESI 2002 April; (O10) ESI 2002 December; (Oil) ESI 2003 February; (012) ESI 2004 March; (013) ESI 2005 January; 
(014) LRIS 1999 August; (015) LRIS 2000 February; (016) LRIS 2003 June; (017) Kast 1998 December; (018) Kast 2000 September; (019) Kast 2001 June; (O20) 
SDSS Spectrum; (021) Double Spectrograph 2006 April 

Near-infrared spectra: (IR1) SpcX 2002 May; (IR2) Spex 2002 October; (IR3) Spcx 2003 April; (IR4) Spex 2003 September; (IR5) Spcx 2004 May; (IR6) Spcx 2004 
June; (IR7) Spcx 2004 September; (IR8) TIFKAM 2002 April 

I 

References. — (B01) Becker et al. (2001); (Gr96) Gregg et al. (1996); (Gr02) Gregg et al. (2002); (G104) Glikman et al. (2004); (Lac02) Lacy et al. (2002); (Law96) ^ 
Lawrence ct al. (1996); (P98)Perlman ct al. (1998); (S02) Smith ct al. (2002); (V96) Vcrmculcn ct al. (1996); (X94) Xu et al. (1994) ° 



Table 2. F2M Candidates with J SNR < 3 and J - K < 1.7 



R.A. 


Dec 




F u 




B 


R 


j 


H 


K s 


J — K 


R — K 


Tvno 


Rodshift 


O Vi ^ p r 1" i o ti ^ 


Ref 


c 120001 


(J2000) 


(mjy) 




( tti air] 


( ma c i 


1 TTl f\ D" 1 


( tti air] 


(ma (t i 




( TTl PI CI 










HI 


(21 




(31 




(51 














C121 


CI 31 




("1 51 


02 15 42.02 


-02 22 


57.1 


359.46 


403.64 


21.83 


20.38 


>16.42 


16.225 


15.059 


>1.36 


5.32 


QSO 


1.178 




D97 


07 55 27.61 


+24 46 


44.5 


1.13 


1.25 


20.33 


19.36 


17.33 a 


15.885 a 


15.285 a 


2.05 


4.08 










08 20 25.97 


+47 28 


03.1 


1.24 


1.29 


21.09 


19.36 


>16.74 


16.141 


15.117 


>1.62 


4.24 










08 23 24.76 


+22 23 


03.3 


1919.32 


2163.90 


21.07 


19.73 


>16.98 


16.408 


15.418 


>1.56 


4.31 


BLLac 


0.951 


O20 


S93 


08 46 22.29 


+18 22 


50.0 


13.21 


13.67 


21.18 b 


19.42 b 


>16.25 


15.986 


14.871 


>1.38 


4.55 










09 03 55.86 


+20 20 


19.0 


2.26 


5.21 


22.50 


19.48 


>16.85 


16.415 


15.361 


>1.49 


4.12 










09 19 46.61 


+22 58 


22.1 


2.27 


2.56 


22.00 


20.24 


>16.63 


15.458 


15.436 


>1.19 


4.80 










09 23 57.10 


+15 38 


29.5 


1.20 


0.72 


>22.50 


19.20 


>16.64 


15.803 


15.022 


>1.62 


4.18 










09 24 22.34 


+32 23 


55.1 


5.93 


5.63 


21.58 


19.48 


>16.47 


16.392 


15.309 


>1.16 


4.17 










09 52 13.18 


+25 00 


59.0 


11.71 


12.61 


21.14 


18.92 


> 16.51 


15.777 


14.895 


>1.61 


4.02 










10 00 06.47 


+45 49 


50.0 


1.99 


1.52 


21.13 b 


19.82 b 


>16.52 


16.111 


15.182 


>1.34 


4.64 










10 14 31.77 


-02 20 


53.4 


18.02 


27.44 


21.08 


19.22 


17.27 a 


16.079 a 


14.983 a 


2.29 


4.24 










10 20 54.07 


+16 27 


02.7 


2.14 


5.46 


22.12 


19.43 


>16.86 


15.897 


15.176 


>1.68 


4.25 










11 06 22.48 


+20 21 


35.6 


2.34 


2.25 


21.62 


19.67 


> 16.73 


15.839 


15.306 


>1.42 


4.36 








h 


11 40 32.29 


+24 01 


18.0 


3.85 


3.69 


21.23 


19.33 


>16.74 


15.996 


15.047 


>1.69 


4.28 










11 42 07.75 


+15 47 


54.1 


270.14 


283.71 


20.06 


19.48 


>16.36 


15.498 


15.305 


>1.06 


4.17 










12 08 27.76 


+17 08 


14.1 


2.26 


5.21 


22.11 


19.40 


>16.62 


15.263 


15.322 


>1.30 


4.08 


Galaxy 


0.549 


015 


G104 


12 29 24.80 


+25 09 


21.1 


2.07 


1.80 


22.50 


20.21 


>16.78 


16.338 


15.434 


>1.35 


4.78 








G104 


12 40 23.71 


+43 53 


18.4 


1.06 


1.66 


20.43 


19.01 


>16.68 


16.191 


15.008 


>1.67 


4.00 










14 21 12.88 


+36 08 


51.0 


1.29 


1.27 


20.09 


19.85 b 


>15.95 


16.446 


15.002 


>0.95 


4.85 










14 34 34.95 


+29 37 


05.1 


2.40 


1.98 


21.78 


19.39 


> 16.81 


16.189 


15.142 


>1.67 


4.25 










16 07 43.42 


+14 30 


37.6 


1.80 


2.18 


21.44 


19.57 


>16.85 


15.562 


15.293 


>1.56 


4.28 











a J, H, and K s magnitudes are from the 2MASS All-Sky data Release catalog. 

h B and R magnitudes were undetected in GSC-II. We replace these with r and g magnitudes from SDSS. 
Note. — Reference codes for the observations are the same as in Table 1 
References. — (D97) Drinkwater ct al. (1997); (G104) Glikman ct al. (2004); (S93) ? 



Table 3. Line Diagnostics of Narrow-Line Spectra 



Object 


Log([0 III1/H3) 


LogffN III /Ha) 


LoffffS III /Ha) 

OVl if / 


Log([0 11 /Ha) 

o\l J / J 


Type 1 


Type 2 


Type 3 


(1) 


(2) 


(3) 


(4) 


(5) 


f6) 


(7) 


(8) 


0701p5230.ms.dat 


-0.01 


-0.63 


-0.35 


-0.88 


starb 


starb 


agn 


0730p2538.ms.dat 


0.73 


-0.68 


-1.02 


-1.91 


agn 


starb 


starb 


0738p2141.ms.dat 


1.31 


-0.29 


-0.61 


-1.18 


agn 


agn 


agn 


0738p3156.ms.dat 


1.04 


-0.41 


-0.86 


-1.53 


agn 


agn 


agn 


0810pl603.ms.dat 


-0.02 


-0.22 


-0.56 


-1.31 


starb 


starb 


starb 


0817pl958.ms.dat 


-1.10 


-1.07 


-1.60 


-1.55 


starb 


starb 


starb 


0818pl305.ms.dat 


1.00 


-0.35 


-0.54 


-1.16 


agn 


agn 


agn 


0849pl439.ms.dat 


-0.07 


-1.62 


-0.61 


-1.47 


starb 


starb 


starb 


0915pl819.ms.dat 


-0.16 


-0.39 


-0.47 


-1.28 


starb 


starb 


starb 


0915p2850.ms.dat 


1.42 


0.23 


-0.08 


-0.83 


agn 


agn 


agn 


0919pl531.ms.dat 


-0.09 


-0.33 


-0.44 


-1.32 


starb 


starb 


starb 


0946pl840.ms.dat 


0.63 


-0.95 


-0.99 


-1.74 


starb 


starb 


starb 


0949p2325.ms.dat 


1.76 


-0.24 


-0.31 


-1.33 


agn 


agn 


agn 


1004p2943.ms.dat 


0.74 


-0.38 


-0.77 


-1.49 


agn 


agn 


agn 


1015pl207.ms.dat 


0.98 


-0.46 


-0.28 


-0.88 


agn 


agn 


agn 


1018p2135.ms.dat 


0.50 


-1.24 


-1.09 


-1.58 


starb 


starb 


starb 


1022pl628.1d.dat 


0.98 








agn a 


agn a 


agn a 


1022pl929.ms.dat 


0.77 


-0.68 


-0.82 


-1.31 


agn 


agn 


agn 


1139p2743.ms.dat 


0.47 


-0.56 


-0.41 


-0.82 


starb 


agn 


agn 


1202p2615.ms.dat 


1.10 








agn a 


agn a 


agn a 


1654p3435.ms.dat 


0.99 


-0.69 


-1.03 


-1.53 


agn 


agn 


agn 


2328mll07.1d.dat 


0.63 


-0.11 


-0.45 


-0.99 


agn 


agn 


agn 



a These objects were at redshifts higher than 0.55 and Ha was outside our spectroscopic range. Nevertheless, we classify these 



sources as AGN since their [O IIIJ/H/3 ratios were above the extreme starburst dividing line shown in Figure 2 



Note. — Type 1 is the classification determined from the [N II] /Ha versus [O IIIJ/H/3 diagnostic. 
Type 2 is the classification determined from the [S II] /Ha versus [O IIIJ/H/3 diagnostic. Type 3 is the classification determined 
from the [O I] /Ha versus [O IIIJ/H/3 diagnostic. 



I 

CO 

I 



Table 4. Extinction Parameters for F2M Quasars 



E(B - V): Balmer Decrement E(B - V) 



R.A. 


Dec 


z 


K s 


A K 


R-K 


Total line 


Broad line 


Cont. Fit 


Spectrum Used 


(J2000) 


(J2000) 




(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


10 04 24.87 


+12 29 22.4 


2.658 


14.54 


1.09 


6.26 






0.48 


Comb. 


15 49 38.80 


+12 45 07.9 


2.373 


13.52 


0.65 


4.09 


0.188+0.012 


0.141+0.013 


0.26 


Single:IR 


15 31 50.48 


+24 23 17.6 


2.287 


14.69 


0.38 


4.48 


-0.005+0.046 


0.170+0.073 


0.16 


SingledR 


22 22 52.78 


-02 02 57.4 


2.252 


15.14 


0.90 


4.90 


-0.632+0.038 




0.39 


Comb. 


01 34 35.68 


-09 31 03.0 


2.220 


13.55 


1.08 


7.25 


0.203+0.005 




0.58 


Comb. 


14 27 44.34 


+37 23 37.4 


2.168 


15.07 


0.60 


4.52 


0.072+0.024 




0.27 


Comb. 


07 38 20.10 


+27 50 45.5 


1.985 


15.26 


0.91 


5.54 






0.50 


Comb. 


09 21 45.69 


+19 18 12.6 


1.791 


14.55 


1.20 


5.57 






0.65 


Opt. 


13 44 08.31 


+28 39 32.0 


1.770 


14.78 


0.96 


6.02 






0.53 


Opt. 


10 36 33.54 


+28 28 21.6 


1.762 


15.27 


0.85 


4.74 






0.47 


Opt. 


13 59 41.18 


+31 57 40.5 


1.724 


14.79 


0.82 


6.01 






0.47 


SingledR 


13 41 08.11 


+33 01 10.3 


1.720 


14.85 


0.98 


5.95 






0.56 


Comb. 


09 06 51.52 


+49 52 36.0 


1.635 


15.14 


0.16 


5.66 






0.10 


Single:Opt. 


09 16 08.57 


-03 19 39.9 


1.560 


14.82 


0.52 


4.72 


0.917+0.075 




0.33 


Comb. 


10 46 42.55 


+16 04 17.4 


1.440 


15.47 


0.60 


4.47 






0.41 


Single:Opt. 


13 53 08.65 


+36 57 51.2 


1.311 


14.26 


1.20 


6.54 






0.90 


Comb. 


09 16 48.91 


+38 54 28.3 


1.267 


14.75 


0.17 


4.02 






0.13 


Single:Opt. 


09 27 44.38 


+39 30 38.4 


1.170 


14.59 


0.77 


4.11 






0.65 


Single:Opt. 


11 46 58.31 


+39 58 34.2 


1.088 


13.53 


-0.04 


4.64 






-0.03 


Single:Opt. 


09 04 50.51 


-01 45 24.5 


1.005 


14.90 


1.10 


5.46 


0.404+0.110 




1.06 


Comb. 


07 29 10.35 


+33 36 34.0 


0.957 


14.52 


0.83 


5.64 






0.84 


Opt. 


10 12 30.49 


+28 25 27.2 


0.937 


15.27 


0.62 


5.53 






0.46 


Comb. 



Table 4 — Continued 



E(B - V): Balmer Decrement E(B - V) 



R.A. 


Dec 


z 


K s 


A K 


R-K 


Total line 


Broad line 


Cont. Fit 


Spectrum Used 


(J2000) 


(J2000) 




(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 




' (1) 


' (2) ' 


(3) 


(4) 


(5) 


(b) 


(7) 


(8) 


(9) 


(10) 


I 1 r\n nT or 

II Ob 07.2b 


i OO 1 O A'l C\ 

+28 12 47.0 


n o a o 

0.842 


1 A P O 

14. b3 


0.17 


tr on 

5.39 






n 1 n 

0.19 


bmgle:Upt. 


nn i r ni 171 

(jy 15 ui.71 


i o A 1 O 1 O O 

+24 18 12.2 


n o a o 

0.842 


1 O r ~7C\ 

13.79 


n o c\ 

O.bO 


7.01 






n r*o 

0.68 


Comb. 


no or no nr 

(Jo 25 02.05 


i A'l 1 / • r* o n 

+47 lb 52.0 


n ono 

0.803 


1/111 

14.11 


O.bl 


O O A 

b.34 






n 71 

0.71 


Comb. 


11 tr 1 O /I nT 

11 51 24. 07 


i ro CT n C7 /I 

+53 59 57.4 


n 'von 

0.780 


i cr i n 

15.10 


C\ O A 

0.b4 


5.70 






n TP 

0.7b 


Comb. 


1 o ot no 01 

12 27 Do. 21 


i r n co o 

+50 53 5b. 3 


n 'ifo 

0.7b8 


1 A f* 1 

14. bl 


n oo 

0.32 


/I 1 /' 

4.1b 


n <vv i n r 

— 0.22b±0.0b5 




n o o 

0.38 


Comb. 


1 r> r f apt -t 1 

lb 5b 47.11 


i OO 01 OP >"7 

+38 21 3b. 7 


n wo 

0.732 


1 rr nn 

15.09 


n /in 

0.49 


5.71 






n a 1 

O.bl 


(Jpt. 


ii r n o i o a 

11 59 31.84 


i on i/i /ion 

+29 14 43.9 


0.730 


11 /IO 

11.48 


n n a 

—0.04 


P 1 1 

b.ll 






n n tr 

—0.05 


bmgle:(Jpt. 


1 a nn o a tr r* 

lb 00 34.56 


1 O tr OO <V7 n 

+35 22 27.0 


n fm 7f\ f ~7 

0.707 


1 A 1 *7 

14.17 


n i o 

0.18 


/I A O 

4.42 






n oo 

0.22 


Single: IK 


11 1 O 1 1 f\f* 

11 18 11. Ub 


nn o o ^i n 

—00 33 41.9 


0.b8b 


1 /I CT O 

14.58 


n o i 

O.bl 


tr on 

5.29 






n 7n 

0.79 


Upt. 


11 io r A r o 

11 13 54. bo 


i 1 O A A o o n 

+ 12 44 38.9 


n ^01 


1 O 

13. bb 


n o o 

0.88 


rr o o 

5.33 






i i tr 

1.15 


(Jpt. 


oo on no oo 

23 39 03.83 


nn i o oi n 

—09 12 21.0 


n rr n 

O.bbO 


i /i nn 

14.09 


n on 

0.80 


/I tr A 

4.54 






1 nr 

1.0b 


Upt. 


1 r a o ro on 

lb 42 58.80 


1 On AO O T o 

+39 48 37.2 


n tr no 

0.593 


1 O ' ) / ' 

12.3b 


n no 

0.02 


A OO 

4.32 






n no 

0.03 


Comb. 


23 25 49.88 


-10 52 43.4 


0.5b4 


15.24 


0.34 


4.46 


0.408+0.005 


0.556+0.010 


0.50 


Comb. 


15 32 33.19 


+24 15 2b.8 


0.5b2 


14.9b 


0.48 


5.13 






0.70 


Comb. 


08 41 04.98 


+3b 04 50.1 


0.553 


14.92 


0.71 


5.88 






1.05 


Comb. 


08 10 58.98 


+41 34 02.b 


0.50b 


14.22 


0.02 


4.11 


-0.159+0.008 


0.245+0.009 


0.04 


Comb. 


08 34 07.01 


+35 0b 01.8 


0.470 


14.65 


0.4b 


4.31 






0.74 


Opt. 


11 17 33.83 


-02 3b 00.3 


0.4b3 


14.05 


0.17 


4.70 


0.121+0.013 


0.032+0.015 


0.29 


Single:Opt. 


lb 18 09.72 


+35 02 08.5 


0.447 


14.05 


0.b5 


4.20 






1.08 


Single:IR 


02 29 50.bb 


-08 42 35.5 


0.440 


14.84 


0.09 


4.34 


0.233+0.015 


0.362+0.028 


0.15 


Opt. 


08 43 32.28 


+49 21 lb. 2 


0.420 


15.13 


0.22 


4.37 


-0.106+0.018 


0.176+0.021 


0.38 


Single:Opt. 


14 15 22.83 


+33 33 Ob.5 


0.41b 


14.30 


0.35 


4.b2 






O.bO 


Single:IR 



Table 4 — Continued 



R.A. Dec z K s A K R-K 

(J2000) (J2000) (mag) (mag) (mag) 

' (1) ' (2) ' (3) (4) (5) (6) 



08 30 


11 


13 


+37 


59 


51 


.9 


0. 


.414 


14. 


.58 


0. 


43 


4. 


22 


17 08 


02 


.53 


+22 


27 


25 


.7 


0. 


.377 


14. 


.57 





.59 


4. 


.77 


17 24 


28. 


.03 


+33 


31 


19. 


.9 


0. 


.370 


14. 


.46 


0. 


.51 


4. 


.02 


12 09 


21 


.17 


-01 


07 


17. 


.0 


0. 


363 


13. 


.72 


0. 


.49 


4. 


14 


14 08 


11 


.61 


+32 


43 


50. 


.2 


0. 


.338 


15. 


.35 


0. 


.54 


4. 


.24 


13 04 


04. 


.07 


+46 


12 


53 


.6 


0. 


.315 


14. 


.32 


0. 


.77 


4. 


00 


17 04 


25. 


.31 


+30 


37 


08. 


.7 


0. 


.311 


14. 


.92 


0. 


.28 


5. 


88 


00 36 


59. 


.85 


-01 


13 


32 


.3 


0. 


.294 


13. 


.65 


0. 


.40 


5. 


.89 


07 53 


41 


.11 


+25 


06 


39 


.4 


0. 


.292 


15. 


.22 


0. 


.51 


4. 


20 


13 07 


00 


.62 


+23 


38 


05 


.3 


0. 


.275 


13. 


.45 


0. 


.28 


5. 


92 


22 16 


33 


.73 


-00 


54 


51, 


.2 


0. 


.200 


13. 


.78 


0. 


.25 


4. 


13 


08 17 


39 


.57 


+43 


54 


20 


.0 


0. 


.186 


14. 


.09 


0. 


.32 


4. 


.57 



E(B - V): Balmer Decrement E(B - V) 



Total line 


Broad line 


Cont. Fit 


Spectrum I 


(mag) 


(mag) 


(mag) 




(7) 


(8) 


(9) 


(10) 


0.650+0.008 


0.661+0.013 


0.74 


bmgle:Upt. 


2.408+0.130 


3.480+0.566 


1.05 


Comb. 






0.91 


Opt. 






0.89 


Single:IR 


1.294+0.035 


1.855+0.156 


1.01 


Single:Opt. 






1.49 


Single:IR 






0.54 


Single:IR 


1.313+0.027 


2.544+0.256 


0.79 


Single:Opt. 






1.02 


Single:Opt. 


1.106+0.031 


1.736+0.101 


0.58 


Single:Opt. 


0.416+0.040 


0.406+0.054 


0.55 


Comb. 






0.73 


Opt. 
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Table 5. QSO Counts and Surface Densities in if-Band 



K Magnitude Range 


FBQS II 


FBQS III 


uvx a 


F2M 


K 




(deg- 2 


0.5 mag -1 ) 





11 


- 


11 


5 


0.0011±0 


0006 






0.0004±0.0004 


11 


5 - 


12 





0.0004±0 


0004 


0.002±0.002 






12 


- 


12 


5 


0.003±0 


001 






0.0004±0.0004 


12 


5 - 


13 





0.004±0 


001 








13 


- 


13 


5 


0.006±0 


002 






0.0004±0.0004 


13 


5 - 


14 





0.014±0 


002 


0.007±0.003 


0.007±0.005 


0.0029±0.0011 


14 


- 


14 


5 


0.030±0 


003 


0.031±0.007 


0.04±0.01 


0.0041±0.0012 


14 


5 - 


15 





0.044±0 


004 


0.042±0.009 


0.07±0.02 


0.007718:885? 


15 


- 


15 


5 


0.052±0 


004 


0.07±0.01 


0.12±0.04 


o.oo52±S:8SS 


15 


5 - 


16 





0.030±0 


003 


0.07±0.01 


0.04±0.02 




16 


- 


16 


5 


0.002±0 


001 


0.010±0.004 









(deg" 2 ) 


(deg" 2 ) 


(deg" 2 ) 


(deg -2 ) 


- 15.5 


0.154±0.008 


0.15±0.02 


0.23±0.05 


o.o2H8;88^ 



a This sample is made up of the LBQS and 2QZ samples combined as described in 
§5.1 of Paper I 
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Table 6. 3.6 cm and 20 cm VLA Observations 



FIRST VLA 20 cm a NVSS VLA 3.6 cm a 



Object 
















ai.4GHz/8.3GHz 






(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 




(1) 




(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


r Z1V1 J U ( Zy 1U. 


OtOOOUOI 


Q OR 
O.zO 


q 17 
o.l / 


9 Q1 

z.ol 




9 K 

z.o 


n kr 
u.oo 


n 

u.oo 


n 70 
— u. < y 


F9A/T 707^890 
r zivu u * oozu. 


. 1-f-Z < ou4o 


9 fi/l 

Z.04 


11 
Z.oo 


1 /in 

1 .4U 




"X 7 


n A7 

U.4 / 


n 'zr 

U.oO 


n R^ 
— u.oi 


r zivijuoiuoo. 


Q_l_/U ^409 
,y-r z ±i04uz 


1 fiO 08 
IDy .yo 


1 80 Q/1 

ioy .o4 


918/19 
Zlo.4z 


01 Q O/l 


910 7 

ziy . i 


OCiK f\R 

zuo.uo 


9no ki 
zuy .01 


n n/i 

— U.U4 


F9A/T Tftfil 7^Q 
r zivuuoi / oy. 


c:_i_/i q c.yi on 


q 1 o 
o.lU 


9 80 
Z.8U 


9 ^9 
Z.OZ 


z.oo 


O 1 
Z. 1 


n a o 
u.4y 


n ro 

U.Oz 


n 09 
— u.yz 


F9A/T 7089^.09 
r ZlVlJUoZOUZ. 


n I /i7i Rf;ob 
.U-f-4 / 100Z 


£1 1 9 
01. lz 


R1 OA 
0O.Z4 


^1 8^ 
01. 80 


KO 71 
OZ. ( 1 


^n k 
ou.o 


1 ^ £8 
10. 0o 


1 K OQ 

10. yo 


n £7 

— U.O ( 


F9A/T in&^ni i 
r ziviJUoouii. 


.i-t-o i oyoi 


RAO 
0.4Z 


R CO 
0.0U 


R 70 

D. ( y 


4.04 


K A 
0.4 


1 09 

l.Uz 


n 07 
u.y / 


1 C\R 
— L.V0 


F9A/T 708^/107 


.UtoOuOUI 


1 99 

l.zz 


n kk 
u.oo 


O 79 

u. ( z 




1 A 
1.4 


n 1 9 

U. lz 




1 nn 

— l.UU 


F9A/T 708/11 0/1 
r zivu uo^iu'i. 


y-rouu4ou 


R A O 


R 

0.08 


R RC\ 


D.oy 


7 A 
1 .4 


Q 09 

o.yz 


A f\R 
4.U0 


n 90 
— u.zy 


F9A/T 7000/1^0 
r zivijuyu40U. 


.0 — U140Z^± 


9 9^ 
z.zo 


9 7^ 
Z. / 


9 n*? 


/i 


9 8 
Z.o 


n aa 

U.44 


n "^o 
u.oy 


n afi 
— U.oO 


F9A/7 TOQOfi^l 
J? zivuuyuuoi. 


.on^^tyozou 


£7 87 


7c; 1 1 
I O.ll 


( . 10 


Ot.iD 


£7 1 
( .1 


78 7R 
lo. 10 


8n nn 

oU.UU 


n no 
u.uy 


F9A/TT0Q 1 ^01 
r zivuuyioui. 


7_l_9zL1 81 9 


q an 
y .ou 


1 n 1 9 

1U.1Z 


K 8/1 
O.o4 


^ fi7 
0.0 * 


o n 
y .u 


oa nn 

Z4.UU 


oa z.r\ 

Z4.0U 


n 70 
u. / y 


F9A/T 70Q1 fiOfi 
r zivijuyiouo. 


.o — uoiyoy 


117 OK 
LL 1 .ZO 


19^ 77 
IzO. 1 1 


1 1 8 /in 
llo.4U 


1 1 Q C/| 

llo.04 


1 97 A 
LZ 1 .4 


KA 97 
04. Z ( 


on 
oo. yu 


n /i a 

— U.44 


F9A/T 70Q91 /L^ 
r zivuuyzi40. 


fill Q1 SI 9 


/i /in 

4.4U 


A 71 
4. / 1 


Q 8/1 


/I 1 O 
4.1y 


A 7 
4. / 


1 OQ 
l.Uo 


i 1 n 

1.1U 


n 7/i 

— U. ( 4 


F9A/T 71 00/19/1 
r zivij 1UU4Z4. 


Q_i_1 99Q99 
,o~rizzyzz 


1 1 AO 
11. 4Z 


1 9 Q9 
Iz.oZ 


1 1 Qfi 
11. oD 


1 1 fiO 
1 1 .0z 


118 
11.8 


I 1 09 

I I .yz 


19 10 

iz.iy 


n nQ 

U.Uo 


F9A/T T1 01 9^0 
r zivu luizou. 


zl 4-989^.97 


O 9Q 

y .zo 


8 7fi 
o. 10 


71 

0. ( 1 


1 1 /in 
1 1 .4U 


inn 

1U.U 


8 QO 
o.oU 


8 KO 
o.Oz 


n 91 

U.zl 


F9A/T 71 0^*3*3 
r ZIVIJ lUoOoo. 


c._i_989891 

.OtZoZOZI 


/I 9p; 
4.Z0 


A Q7 
4.o ( 


Q 07 

o.y ( 


o.O * 


C 7 
0. / 


/I fi8 
4.08 


A 7K 


n no 
u.uy 


F9A/T T1 1 0fi07 


9 i 981 9/17 

ZTZOlZ'i / 


911 A 7 
Zl 1 .4 f 


91 A 70 
Zl4. ( y 


99/1 0/1 
ZZ4.y4 


oofi ofi 
zzO.zO 


99/1 7 
ZZ4. 1 


Qan fifi 
ooU.OO 


Q8fi 98 
oSO.Zo 


n 90 
u.zy 


F9A/TT1 1 1 QC/i 
r zivu 1 11004. 


9zLzL^8 b 


9 OO 

z.yy 


A 9Q 
4.ZO 






A R 
4.0 


9 no 
z.uy 


9 0^ 
Z.UO 




F9A/T T1 1 1 7^*3 
r ZIVIJ 111 * oo 


s n9 c ffinnt> 

o — UZoOUU 


fine n/i 

0U8.U4 


Ofi1 OR 

yoi .zo 


on^ 
yuo.oo 


yoU. lo 


OOfi 1 

yyo.i 


9QQ 00 

zoo. yy 


979 A 8 
z / z.4o 


n 7£ 

— U. ( 


T?01\/T T1 1 1 01 1 
r ZIVIJ llloll. 


n nnQQ/11 
U— UUoo41 


1 qo 

l.oU 


1 81 
1.81 


n 78 

U. ( o 




9 8 
Z.o 


n 98 

U.ZO 


n 99 
U.zz 


n ^7 

— U.O ( 


F2MJ115124 


,0+535957 c 


3.5z 


3.60 


o on 
3.29 




3.5 


0.55 


0.71 


1 AA 

— 1.00 


F2MJ115931. 


,8+291443 


1855.80 


1952.59 


2023.5 


2027.4 


2030.8 


1743.2 


1766.3 


-0.08 


F2MJ120921. 


.2-010717 


1.42 


1.25 


1.77 




2.8 


1.11 




-0.26 


F2MJ122703. 


.2+505356 


3.87 


4.27 


3.66 


3.94 


3.9 


0.99 


0.94 


-0.73 


F2MJ130404. 


.0+461253 


2.14 


2.63 


2.38 


2.59 


3.2 


.75 




-0.65 


F2MJ130700. 


,6+233805 


2.99 


2.63 


2.68 


2.38 


3.0 


0.87 


0.79 


-0.63 


F2MJ134108. 


.1+330110 


69.41 


69.81 


63.12 


63.28 


71.6 


23.06 


23.50 


-0.56 


F2MJ134408. 


,3+283932 


9.76 


10.41 


7.95 


10.41 


10.8 


1.82 


1.80 


-0.83 


F2MJ135308. 


,6+365751 


3.71 


3.32 


3.75 


3.24 


3.5 


0.76 


0.69 


-0.89 


F2MJ140811. 


,6+324350 


1.09 


2.26 


0.63 




1.4 


0.69 




0.05 


F2MJ141522. 


,8+333306 b 


5.21 


6.21 


4.12 


5.51 


8.9 


1.32 




-0.64 


F2MJ142744. 


,3+372337 


1.71 


1.52 


0.99 




1.4 


0.96 




-0.02 


F2MJ153150. 


,4+242317 


2.24 


2.00 


6.66 




1.4 


0.28 




-1.78 


F2MJ153233 


l+241526 b ' c 


7.43 


9.37 


11.69 




5.3 


7.39 


7.54 


-0.26 
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Table 6 — Continued 



FIRST VLA 20 cm a NVSS VLA 3.6 cm a 



Object 




F int 










f \nt 


a 1.4GHz/8.3GHz 




(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


F2MJ154938. 8+124507 


1.80 


1.79 


1.94 


2.07 


2.0 


0.53 


0.49 


-0.73 


F2MJ160034.5+352227 


3.17 


3.30 


1.17 




3.4 


0.60 


0.67 


-0.37 


F2MJ161809. 7+350208 


13.73 


16.22 


34.25 


37.32 


35.8 


7.02 


8.50 


-0.89 


F2MJ165647.1+382136 


4.12 


3.89 


4.09 


4.34 


4.5 


0.71 


0.72 


-0.98 


F2MJ170802. 5+222725 


1.31 


1.74 


1.05 




2.8 


0.37 


0.37 


-0.58 


F2MJ172428.0+333119 b 


1.58 


1.35 


<1.20 


<1.20 


2.8 








F2MJ221633. 7-005451 


1.31 


0.80 


1.44 




2.0 


0.21 


0.32 


-1.08 


F2MJ222252. 7-020257 


2.43 


2.66 






2.5 


0.16 


0.30 




F2MJ232549. 8- 105243 


1.46 


1.10 






2.8 


0.25 


0.20 




F2MJ233903. 8-091221 


4.34 


4.06 


2.70 




3.9 


0.35 


0.39 


-1.14 



Observations were conducted using the CnD configuration of the VLA, unless otherwise noted. 
Observations were conducted using the D configuration of the VLA 

c Gaussian fitting was forced to a point source, with the FWHM fixed at the synthesized beam size. 



Table 7. Optical Properties of Blue F2M Quasars 











GSC II 






APM 






SDSS 








RA 


Dec 


z 


F 


J 


J - F 


E 


O 


O-E 


r 


9 


g-r 


K 


E(B-V) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


08 10 58.98 


+41 34 02.6 


0.506 


18.33 


18.27 


-0.06 


18.27 


18.67 


0.40 


18.64 


18.53 


-0.11 


14.199 


0.06 


09 06 51.52 


+49 52 36.0 


1.635 


>20.80 


20.82 


<0.02 


>20.00 


>21.50 




21.23 


21.50 


0.27 


15.141 


0.10 


11 46 58.31 


+39 58 34.2 


1.088 


18.17 


18.47 


0.30 


18.03 


19.02 


0.99 


18.65 


18.30 


0.35 


13.530 


-0.04 


11 59 31.84 


+29 14 43.9 


0.730 


17.59 


17.67 


0.08 


16.58 


17.53 


0.95 








11.474 


-0.06 


16 42 58.80 


+39 48 37.2 


0.593 


16.68 


17.42 


0.74 


14.70 


15.48 


0.78 


16.47 


16.57 


0.10 


12.398 


0.03 



Table 8. Radio Properties of Blue F2M Quasars 



FIRST Spectral Index 



RA 


Dec 






NVSS 


6 cm 


n FIRST 
"6cm 


n NVSS 
"6cm 


logK a 


NED 


Comment 


(J2000) 


(J2000) 


(mjy) 


(mjy) 


(mjy) 


(mjy) 








(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


08 10 58.98 


+41 34 02.6 


170.0 


189.3 


219.7 


168 


-0.10 


-0.22 


2.98 


y 


RASS X-ray source Brinkmann et al. (2000) 


09 06 51.52 


+49 52 36.0 


67.9 


75.1 


57.1 


96 


0.20 


0.42 


3.51 


y 




11 46 58.31 


+39 58 34.2 


335.6 


338.5 


330.9 


733 


0.62 


0.64 


3.55 


y 




11 59 31.84 


+29 14 43.9 


1855.8 


1952.6 


2030.8 


1542 


-0.19 


-0.22 


3.64 


y 


Blazar Ghosh et al. (2000) 


16 42 58.80 


+39 48 37.2 


6050.1 


6598.6 


7098.6 


6560 


-0.01 


-0.06 


4.20 


y 


Blazar Webb et al. (1994); Raitcri et al. (1998) 



a Radio loudness parameter, defined as in Stocke et al. (1992) logR = log(/(5GHz)//(B)) 
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Fig. 1. — Schematic diagram of our selection process. Major catalogs are indicated by 
dotted-lined boxes and the matching radius is shown in the circle at the node of each match. 
Sources in the final catalog are indicated by bold-lined boxes. 
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Fig. 2. — Line diagnostics for our narrow-line spectra, plotted on diagrams from Kewley 
et al. (2006). The left panel shows the [N II] /Ha versus [O IIl]/H/3 diagnostic. In this panel, 
objects below and to the left of the dotted line are pure starbursts. The middle panel shows 
the [S II] /Ha versus [O IIl]/H/9 diagnostic. The right panel shows the [O I] /Ha versus 
[O IIl]/H/5 diagnostic. Objects to the left and below the solid line are starburst dominated 
(filled cirlcles), while objects above and to the right of the line are AGN dominated (open 
circles). 
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Fig. 3. — All available spectra of F2M candidates identified as quasars, ordered by redshift. 
The dotted lines show expected positions of prominent emission lines in the optical and near- 
infrared: Lya 1216, N V 1240, Si IV 1400, C IV 1550, C III] 1909, Mg II 2800, [O II] 3727, 
R5 4102, H 7 4341, H/3 4862, [O III] 4959, [O III] 5007, Ha 6563, Pa 7 10941, Pa/5 12822. 
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Fig. 3b. — Continued. Spectra of F2M quasars. 
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Fig. 3c. — Continued. Spectra of F2M quasars. 
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Fig. 3d. — Continued. Spectra of F2M quasars. 
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Fig. 3e. — Continued. Spectra of F2M quasars. 
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Fig. 3f. — Continued. Spectra of F2M quasars. 
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Fig. 3g. — Continued. Spectra of F2M quasars. 
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Fig. 3h. — Continued. Spectra of F2M quasars. 
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Fig. 3i. — Continued. Spectra of F2M quasars. 
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Fig. 3j. — Continued. Spectra of F2M quasars. 
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Fig. 4. — Histogram of the K magnitudes for the F2M candidates {unshaded). All spec- 
troscopically identified objects are overplotted in the shaded histogram, and the quasars are 
overplotted in black. 
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Fig. 5. — B — K distribution of F2M quasars compared with FBQS II and III; arrows denote 
bins which are lower limits. 
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Fig. 6. — Redshift distribution of as a function of color of F2M quasars (triangles) compared 
with FBQS II (open circles) and III (filled circles). 
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Fig. 7. — Spatial density of quasars on the sky of F2M red quasars compared with FBQS 
II and III as well as radio-detected UVX-selected quasars from the LBQS and 2QZ quasar 
surveys. 
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Fig. 8. — Examples of a reddened composite spectrum fitted to two spectra. Top - 
F2M003659. 8— 011332 has an optical spectrum and is best fit by a composite spectrum red- 
dened with E(B - V) = 0.79, which is overplotted in red. Bottom - F2M154938.8+124507 
has only a near-infrared spectrum and is best fit by the same composite reddened with 
E(B — V) = 0.26, overplotted in red. 
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Fig. 9. — Two examples of combined optical and near- infrared spectra with the best-fit 
reddened composite overplotted in red. 
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Fig. 10. — An example of a spectrum whose combined optical plus near-infrared spectrum 
was poorly fit by a reddened composite, overplotted in red. The optical spectrum alone 
provides an excellent fit to the reddened composite, overplotted in blue. 
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Fig. 11. — Examples of Ha (left) and H/3 (right) fits for computing Balmer decrements using 
the different modeling techniques described in the text. The object's name is labeled on the 
left panel; the redshift and E(B — V) value is given above the right panel for each object. 
The total fit is plotted with a solid line. When applicable, the narrow-line components 
are overplotted with a dotted line and the broad-line component is marked with a dashed 
line. The linear continuum is overplotted with long-dashed line, (a) High signal-to-noise 
near-infrared spectrum with no narrow-line model, (b) Lower signal-to-noise near-infrared 
spectrum with no narrow-line model, (c) High-signal-to-noise optical spectrum with [O III] 
used for the narrow-line model, (d) High signal-to-noise spectrum with [S II] used for the 
narrow-line model. 
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Fig. 12. — Dereddened .fT-band absolute magnitude as a function of redshift. The colors of 
the circles correspond to the amount of extinction, ranging from no extinction represented 
by light shading (bright yellow) to heavily reddened (dark red). The dotted lines indicate 
the survey limit (K < 15.5) for increasing amounts of extinction. The small dots are FBQSII 
and FBQSIII quasars, which we assume are largely unabsorbed. 
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Fig. 13. — Color excess derived from Balmer Decrement measurements versus color excess 
derived from reddened composite fits to spectra. Square points are from Balmer decrements 
measured from the total line profile, while Balmer decrements derived from the broad com- 
ponent of the lines are plotted with stars. The Balmer decrement in the total hydrogen line 
profile versus the reddening in the continuum is best fit by the solid line. The dashed line 
shows the best fit to the extinctions derived from the Balmer decrement in only the broad 
component of the hydrogen line. The dotted line indicates equivalent measurements between 
the color-excess derived from the Balmer decrement versus the reddening in the continuum. 
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Fig. 14. — Color excess derived from composite fitting (open circles) and from Balmer 
decrements using the total line profile (closed circles) as a function of redshift. Fits from the 
same object are connected with a line. The error bars on the open circles indicate the range 
of values obtained from the steep and flat template spectra (spectral slopes a u = —0.25 and 
a u = —0.76, respectively). 
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Fig. 15. — Extinction curves from Gordon & Clayton (1998) with varying R v values spanning 
the Infrared-though-UV wavelength range. The extinction as a function of wavelength is 
identical in the UV and optical for the different dust laws. They diverge around 7000A 
where larger dust grains characterized by larger R v values lead to less extinction in the 
infrared. 
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Fig. 16. — Color excess as a function of redshift for R v = 3.1 (filled circles). The error 
bars represent the extreme extinctions derived from the Ry = 2.2 and 5.5 dust laws. Higher 
extinction is derived from the Ry = 2.2 dust law, since it is more sensitive in the near- 
infrared, while the lower value is computed with the Ry = 5.5 dust law as it is less sensitive 
to dust at these wavelengths (see Figure 15). 
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Fig. 17. — Flux-density as a function of wavelength for power-law spectra with flat spectral 
indices (F u oc v a with a = —0.5,-0.2,0.0,0.2). The left-hand y-axis presents the 2.2 /im 
brightness in magnitudes, while the right-hand y-axis presents the flux density in milliJansky. 
This model makes the simple assumption that there is no change in the spectral shape over 
six dex in wavelength between 20 cm and 2 /im and determines the amount of flux contributed 
in the if-band if the 1.4 GHz flux density is 1 mJy 
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Spectral Index Spectral Index 

Fig. 18. — Spectral index histogram for 41 F2M and 214 FBQS quasars, normalized by the 
number of quasars in each sample. Left - 11 F2M and 193 FBQS sources brighter than 10 
mJy at 20 cm. Right - 30 F2M and 21 FBQS sources fainter than 10 mJy at 20 cm. 
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Fig. 19. — Spatial density of quasars on the sky of F2M red quasars corrected for K-hand 
absorption and compared with FBQS II and III (assuming no absorption for those quasars). 
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Fig. 20. — Blue quasars detected in the F2M sample, displayed as in Figure 3. The diamond 
symbols represent three SDSS bands, g, r, and i, followed by the three 2MASS bands, J, H, 
and K s . 



